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ABSTRACT 


Long-wave radiation loss maps, based on Explorer VII measurements of terrestrial radiation at night, are analyzed 


and compared with composite nephanalyses and frontal analyses. 


Results indicate a definite relationship between the 


radiation centers and their corresponding surface low and high pressure centers, their locations, 24-hour intensifica- 


tions and movements, and the conformity of these movements to the 500-mb. geostrophic flow. 


applications to analysis and forecasting are noted. 


1. INTRODUCTION 


The radiation-balance experiment aboard the Interna- 
tional Geophysical Year satellite, 1959 Iota (Explorer 
VII) was developed by a University of Wisconsin team 
and is described in [8]. 

The primary purpose of the experiment is to measure 
the solar, reflected, and terrestrial radiation currents to 
obtain the radiant heat flow to and from the earth, and 
ultimately to obtain a clearer understanding of the 
driving force behind the circulation of the atmosphere. 
Observations of the long-wave radiation currents from 
the shadow zone portion of the earth measured by the 
omnidirectional bolometer were used to draw composite 


radiation loss maps. The relationship of these maps to 


other meteorological analyses and parameters and their 


application to analysis and forecasting are described in 
this paper. 
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Some of the potential 


2. LONG-WAVE RADIATION LOSS MAPS 
RADIATION OBSERVATIONS 


The satellite’s orbit is inclined 50.3° to the earth’s 
equator and varies in altitude from 1,100 km. at apogee 
to 550 km. at perigee. During the period studied, early 
December 1959, the satellite came out of the sun at 50° 
N., traveled southeastward and entered the sunlight zone 
in the Southern Hemisphere. The orbital period is 101 
minutes and radiation observations are taken every 42 
seconds, during which the satellite travels 2.5° of latitude 
Limitations of this orbit are the restriction of 


N. the 


of local time along the orbit, and scarcity of observations 


(fig. 1). 
observations to between 50° and 3., variation 
at lower latitudes. 

The spherical sensor integrates the radiation received 
from the entire field of view. The amount received from 
a unit area of the earth is weighted according to its posi- 
tion with respect to the satellite, being inversely propor- 
tional to the square of the distance from the point on 
earth to the satellite and directly proportional to the 
cosine of the angle between the local zenith and a line 
to the satellite as observed from that point [8]. For 
example, at satellite altitude of 700 km. the total radius 

119 
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Figure |.— Long-wave radiation loss map, December 1, 1959. Lines are isolangleys, connecting equal long-wave radiation loss values i 
Langleys per minute * 10-3. Surface pressure centers and fronts are composite to the satellite passage time. Surface pressure centers 
marked by circles, are designated L for Low and H for High, plus an Arabie number (thus L-1). Radiation centers, marked by squares 


are distinguished from surface centers by the inclusion of R (thus L-R-1 is radiation Low number | which is associated with surface 
Low L-1). Satellite observation points and orbit path are shown by dots. 
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Fiqure 2.—Composite nephanalysis, December 1, 1959, based on 0000, 0600, and 1200 Gm surface maps. Overcast (OVC) areas # 
enclosed within solid lines. Broken lines separate the clear to scattered (SCTD) areas from the broken, 5/8 to 7/8 (BRKN). Wai 
lines indicate the actual cloud area edge is unknown due to lack of data. Arrows show regions of maximum winds, from the OQO0O0 6M 
300-mb. chart. 
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of the circular area viewed on the earth is 25.8° of latitude 
(2,855 km.) measured from the subsatellite point. How- 
ever, 50 percent of the radiation comes from a relatively 
small area whose radius is 6.6° of latitude and which 
comprises only 7.6 percent of the total area viewed. 
Nine-tenths of the radiation comes from only 33.2 per- 
cent of the area viewed [4]. As a result, the changes in 
radiation are heavily weighted toward the area just 
below the satellite. 

The surface of the earth or of a cloud radiates essentially 
as a black body. <A portion of this radiation is trans- 
mitted through the atmosphere and the remainder is 
absorbed by atmospheric gases, principally CO,, H.O, 
and Os. In 
fact, most of the energy leaving the top of the atmosphere 
Nocturnal radiation 


The atmosphere itself reradiates energy. 


comes from the atmosphere itself. 
measurements were used in this study because radio 
reception from the satellite is much better at night when 
This results 
In addition, 


there is greater ionospheric transparency. 
in more data and better global coverage. 
calculations for the nighttime ease are much simpler 
because the white hemisphere is measuring only long-wave 
radiation, 

The bolometer is a thin hemispherical shell thermally 
isolated from a mirror large enough so that no part of the 
satellite itself is seen by the bolometer. Since the satellite 
is spinning rapidly on an axis parallel to the mirror surface, 
the hemispherical bolometer acts like a sphere in space. 
The heat transfer by conduction and long-wave radiation 
to the mirror is small but not negligible compared to 
transfer from the earth The 
energy balance equation for a spinning mirror-backed 


radiative and to space. 


MONTHLY WEATHER REVIEW 


A NOt ae 


2e 
AW)7 TUE DEC] 0c- 


500-mb. map, 1200 emt, December 1, 1959. 


hemisphere in space in the shadow zone of the earth 


reduces to 


ap im =2reo Tt, — H = —K(T,,—T,,) (1) 
ous peer SE Ee 1 - 


8 a a Of a 


where ? ¢ is the total radiation received by the bolometer, 
8 is the solid angle to the earth, @ and ¢ are the infrared 
absorptivity and emissivity of the sensor, o the Stefan- 
Boltzmann constant, H a term equivalent to the time 
constant of the sensor, and A the thermal conductivity 
between the temperatures of the sensor (7.3 and the 
mirror (T,,). 

Actually H7 is a function of temperature because the 
heat transfer process, principally radiation, is temperature 
dependent. KA is also temperature dependent because 
there is radiation exchange between the inside of the 
hemisphere and the mount. However, in order to sim- 
plify the hand calculations, both /7 and K are assumed to 
be constant. Even so, since a=e from Kirchoff’s law, 
the value of the first term in the brackets, the greatest of 
the three, is determined by the absolute temperature 7). 
A most pessimistic value of its error would be +1 °C 
a more likely value is about one-fourth of that. However, 
even an error of 1° C. would give rise to an error of only 
2.5 percent in the final R’ The satellite 
altitude and therefore the solid angle 8 is known to good 
An error of 100 percent in K will give 
Finally, while 


value of 
accuracy. an error 
of 10 percent in the final result. 
librium the second term goes to zero, it might be significant 


at equi- 
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2, 195%. (See fig. 1 for legend 


























Figure 5.—Composite nephanalysis, December 2, 1959. (See fig. 2 for legend. 


just as the satellite comes out of the sun. One is tempted 
to say, therefore, that the radiation values at 50° N. are 
probably too high. Our feeling is that it is best to await 
the more complete calculations being done on a computer. 
However, it will be a surprise if they differ by more than 
10 percent from what we are showing here. In any case 


these results are sufficient to show changes in the patterns 
of outgoing radiation. 

Thus the radiation observation values used in the 
long-wave radiation maps are probably within 10 percent 
of the absolute values. However, the values for three 
maps which do not cover the complete hemisphere should 
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not be used for comparison with average hemispheric 
values. From a synoptic point of view, the radiation 
observations are extremely useful because they are com- 
parable to each other relatively so that the radiation 
patterns and their continuity have meaning. 


ANALYSIS PROCEDURE 


Because the spherical sensor integrates the radiation 
over a large segment of the earth, it is possible to draw 
a composite radiation map using data from successive 
satellite orbit paths, even though there is a time difference 
of | hour and 40 minutes between orbit paths. Lines 
of equal long-wave radiation loss, in Langlevs per minute 

10-* (isolangleys) are drawn on the radiation map, 
using the same principles as those for drawing isobars 
onasurface map or contours on an upper air map. 

Since the radiation maps are not synoptic, composite 
surface frontal and nephanalysis maps, which coincide 
as closely as possible with the times of the radiation 
observations, are constructed based on National Weather 
Analysis Center analyses. The composite nephanalysis 
is drawn separating overcast, broken, and clear to scattered 
considered 


Fog and skv-obscured are 


The lack of data over oceanic and tropical 


areas reports 
as overcast. 
areas is a limiting factor to both composite analyses. 

First the radiation maps were drawn as completely 
and independently as possible. The composite nepha- 
nalvses were then drawn. In the areas north of 35° N. 
to 40° N., where there were sufficient data for independent 
radiation and nephanalyses, the close relationship be- 
Hence the nepha- 


tween the two analyses was apparent. 


nalysis, where available, was used to supplement the 
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500-mb. map, 1200 emt, December 2, 1959. 


radiation analysis in those areas where data were sparse. 
While not best in 
these few, early radiation maps to permit as much study 
as possible of their potential With 
density of radiation observations, a completely independent 


an ideal method, it was considered 


value. increased 
radiation map could be drawn. 

It is interesting to note that in those large sections of 
the oceans and Tropics where there are no surface obser- 
vations for a nephanalysis, the major source of data was 
the satellite radiation observations. Because of the large 
gaps between orbit paths in the Tropics, interpolation 
in drawing isolangleys is difficult. Here 24-hour con- 
tinuity in radiation observations proves invaluable. 


3. RADIATION MAP FEATURES AND RELATIONSHIPS 


The fact that clouds dominate in determining the radia- 
tion pattern is apparent when comparing the radiation 
analyses (figs. 1, 4, and 7) with the nephanalyses (figs. 2, 
5, 8) 40° N. 
analysis. 


in those areas north of where data were 


for their independent Since 


sufficient most 
mid-latitude extensive cloud systems of variable vertical 
development are caused by surface low pressure centers 
and fronts, the resulting radiation map patterns also 
reflect Where the surface Low or front 
is most intense, the cloud tops are higher and colder, 
the 


cloud tops lower, the radiating temperature rises and the 


these sources. 
and the outgoing long-wave radiation is less. As 


amount of radiation measured by the satellite sensor 
increases. On the radiation map this results in a radia- 
tion Low over the area of maximum cloud development, 
and higher isolangley values over the area of lower cloud 
The highest 


tops. radiation readings are found over 








MONTHLY WEATHER REVIEW 


NOVEMBER 1{)}] 

















FIGURE 7. 


Long-wave radiation loss map, December 3, 1959. 


(See fig. 1 for legend.) 




















FIGuRE 8. 


the central and eastern sections of surface Highs where 
clear to scattered conditions generally prevail. 

From the above relationships, one would expect to 
find the radiation low center east of the mid-latitude 
surface low center and the radiation high center east of 
the mid-latitude surface high center; and one does find 


Composite nephanalysis, December 3, 1959. 


(See fig. 2 for legend.) 


this occurring in 24 of the 25 cases for which data are 
available (fig. 11). 

An example of the usual relationship is radiation Low 
L-R-1, located southeast of surface Low L-1 and the 
associated Pacific occlusion on December 3, 1959 (fig. 7). 
It is interesting to note that the orientation of the radiation 
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Figure 10.—500-mb. map, 0000 Gat, December 4, 1959. 
Low and trough parallels the occlusion and area of change in these criteria, 24-hour frontal changes. An 


maximum intensity, and the radiation pattern’s spiral-like 
extension northwestward resembles TIROS cloud photo- 
graphs of occlusions. 


example is the cold front which starts in northwestern 
Canada from the remnants of an old occlusion and extends 
southwestward into the Pacific (fig. 1). The radiation 
map correctly indicates that the front is stronger in the 
region where its clouds produce L-R-4 and the sharp 
radiation trough to its southwest, weaker through the area 


In addition to frontal location, the radiation trough, 
ts sharpness and strength of isolangley gradient, 
indicates frontal intensity, and by the interdiurnal 


by 
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where radiation High H-R-3 reflects only clear to seat- 
tered clouds (fig. 2), and again more intense as it becomes 
a warm front leading into the unstable wave of surface 
Low L-1 where the radiation trough strengthens toward 
L-R-1. The weakening of the radiation trough off the 
Pacific coast on December 2 (fig. 4) reflects the decreased 
activity of the surface cold front. 

The entire radiation trough pattern from L-R-4 to 
L.-R-1 and then southwestward resembles a cyclone 
“family” and, in fact, is the radiation map reflection of one. 
The development of family ‘member’? L-1, from an un- 
stable wave to an occlusion can also be followed on the 
radiation maps during the succeeding days. 

The arrows representing the pattern and location of 
maximum winds on the U.S. Weather Bureau 0000 cour 
300-mb. maps have been plotted on the nephanalyses. 
On December 1, 1959, (fig. 2) the maximum upper wind 
and frontal patterns conformed to that given by [10] and 
[2] for the relationship of the jet stream to a cyclone 
family. Again, the patterns of Lows and 
troughs show the same relationship to the jet stream as do 


radiation 


the surface Lows and fronts of the cyclone family. 

While a radiation Low is usually the reflection of a 
surface Low or front, it can also be caused by a cold or 
snow surface or an extended sheet of high clouds. 
these sources may generally be distinguishable by radiation 
There are no snow surfaces 
indicated in this series of maps. However, it is very likely 
that extended sheets of high clouds caused the tropical 
related to identifiable 


map analysis and continuity. 


radiation Lows which are not 
Such large areas of high clouds in the 
It is felt that if 


above an 


synoptic features. 
Tropics have been reported by [1, 7, 5, 9]. 
there were sufficient radiation observations 
easterly wave, the resulting radiation pattern would show 
a minimum over the high clouds in the convergence area. 
Synoptic data available do not indicate any easterly 


waves for these dates. 
4. 24-HOUR CHANGES IN INTENSITY 


During the intensification of a surface or upper air 
Low there is an increase in size, vertical development, 
and meridional expansion of the cloud pattern. These 
changes are reflected in the radiation pattern. 

The changes in value of the radiation center reflect the 
interdiurnal change in intensity of its surface counterpart 
for five of eight high centers and four of six low centers 
studied here (figs. lla and b). Reasons for the “‘misses”’ 
are difficult to determine but some possiblities are: the 
pressure changes may have been too small to result in 
changes in clouds and radiation intensity; the clouds may 
already have been at maximum height. From the limited 
examples available, one can sav that in all four cases 
when the 24-hour change in magnitude of the surface Low 
was greater than 10 mb., the change in radiation Low 
intensity was in the same sense. 

Intensification of a surface center is also indicated by 


24-hour shifts in the relative positions of its related radia- 
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TABLE 1.—Surface and radiation center motions 


500-mb. geostrophie wind 





24-hr. 24-hr. 
Radiation motion Surface motion 
center (kt.) center (kt.) At radiation | Padiation cen. 
center location ter motion/| ; 
kt (percent 
December 1-2, 1959 
H-R-1 E ll None 
H-R-2 E 22 H-2 ENE 25 WwW 40 55 
H-R-3 SE 19 
H-R-4 NE 20 Hi-4 to H-2 E7 Light 
H-R-5 8 16 H-5 SSE 16 NNW 35 it} 
H-R-6 
H-R-7 ESE 21 H-7 ESE 22 NW 45 45 
L-R-1 NE 2¢ L-1 ENE 24 WSW 45 53 
L-R-2 E 35 L-2 E 35 WSW 60 5s 
L-R-3 E 15 
L-R-4 ENE 26 L-4 ENE 326 SW 65 40 
L-R-5 NNE 19 
L-R-6 E 25 
L-R-7 ESE 15 L-7 SE 14.. W 35 43 
L-R-& 
Decerm ter 2-3, 1959 
H-R-1 
H-R-2 H-2 ENE 32 WNW 50 58 
H-R-3 
H-R-4 H-4 to H-2 NE 22 NW 20 100 
H-R-5 H-5 SE 11 NW 30 33 
H-R-# 
H-R-7 
L-R-1 NE 36 L-1 NE 42 SW 80 45 
L-R-2 
L-R-3 E17 
L-R-4 ENE 31 L-4 E 3. WSW 75 41 
L-R-5 ; 
L-R+ 
L-R-7 7 SE 42 Upper Low extends to 300 mb ’ 
L-R-8 L-&8 to L-7 


Both of 


southeast of surface Low L-1 on December 1, 


3 days of maps. 








tion center. For example, radiation Low L-~R-7 moved 
from southeast to north-northeast of its intensifving sur- 
Limited 


face Low from December 2 to 3, 1959 (fig. 11a). 
data in the Pacific indicate that L—R—1 was located east- 


moved to 





its east-northeast on the 2d, and then to its southeast 
after occlusion took place. On both dates L-1 deepened. 
Reasons for the last change are not readily apparent 
With an intensifying system such as L-1 moving north- 
eastward, one would normally expect some northward 
motion of the radiation Low relative to the surface Low, 
reflecting the movement of highest clouds. However, in 
this case, the radiation Low and trough were apparently 





in the region of most intense convergence and _ vertical 
motion just ahead of the surface position of the occlusion 
(figs. 7, 9, 10). It was also noted in the few available 
cases that radiation Lows related to surface Lows whieh 
were weakening or not changing retained their same rela- 
tive position. L-2 and L-4 and their radiation Lows are 
good examples. 

Another indication of surface Low intensification, based 
on radiation map data alone, is the change in orientation 
of a radiation low pattern from east-west to north-south 


coupled with a decrease in magnitude of the radiation 


This occurred with L-R-1 and L—R-7 over the 
The reverse applies to a weakening sys- 


center. 


tem, as exemplified by L—-R-2 and L-R-4. 


5. 24-HOUR MOTION AND 500-MB. GEOSTROPHIC 


FLOW 


Figures 11a and b and table 1 show a possibly useful 
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Ficure 11.—(a) Surface and radiation center continuity, December 1-3, 1959. The surface pressure center locations are indicated by 


circles, their values in millibars are underlined, and their 24-hour locations are connected by a dashed line. 


tions are indicated by squares, their values are in Langleys per minute 


The radiation center loca- 
10-3, and their 24-hour locations are connected by a solid 


Ine (b 


relationship between the interdiurnal movements of the 
mid-latitude surface and associated radiation Lows and 
Highs. 
Highs available for study was within about 20° and 10 
The seventh 


The 24-hour motion of six of the seven radiation 


kt. of that of their related surface system. 
example, comparing the motions of H-R-4 and H-4, is 


613342 


61 2 


Additional surface and radiation center continuity, December 1-3, 1959. 


complicated by H—4’s becoming absorbed in surface High 
H-2. 


followed that of their surface Low within these same limits. 


The 24-hour movements of all seven radiation Lows 


A positive relationship is found between the 24-hour 
motion of the radiation low center and the average 500- 
mb. geostrophic wind velocity over its area of motion, 
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Computations were made from contour spacings on the 
National Weather Analysis Center 500-mb. maps for 
December | through 4, 1959 (figs. 3, 6,9, 10). The motions 
of four of seven radiation Highs were within 20° of the 
direction and between 40 and 60 percent of the speed of 
the related 500-mb. geostrophic wind velocity (table 1). 

The radiation Lows’ motions show even better correla- 
tion. Sir of the seven cases for which data are available 
were within these limits. The seventh case involves the 
intense surface Low L-7 and radiation Low L—R-7 on 
December 3. A closed Low also developed in the upper 
air maps to at least 300 mb. (highest level available for 
study) and hence the 500-mb. level cannot be considered 
a “steering” map for that system. 

These results indicate that the radiation low center’s 
motion could be used to determine the 500-mb. flow 
pattern, and the approximate location of the corresponding 
surface Low center, its movement, and possible changes in 
its intensity. 

The cirrus clouds causing the tropical radiation Lows 
were probably between 12 and 16 km. Although the 500- 
and 300-mb. analyses extend only to 15°-20° N., the 
contour pattern indicates the westerly flow continued 
farther south at these and higher levels, causing the 
generally eastward motion of H-R-1, L-R-5, and L-R-6. 
L—-R-3 could have been under the influence of the equa- 
torial westerlies which, according to [6], can extend to 
16 km. 
continuity (fig. Ila). 


6. CONCLUSIONS 


The relatively simple omnidirectional white radiation 


These subtropical centers also showed good 


sensor aboard Explorer VII measures long-wave radiation 
at night which, when converted to Langleys per minute, 
can be used to produce radiation maps. In mid-latitudes, 
where adequate data exist, the radiation pattern is consis- 
tent with the nephanalysis pattern and with synoptic 
features. Analysis similar to 
presently used for surface and upper air maps. 


techniques are those 

The radiation map appears to offer many aids to analysis 
and forecasting, particularly for areas of sparse data, 
which should be confirmed by further study. Best 
results are for the mid-latitude low center relationships. 
The radiation low center gives an approximation of the 
location, motion, and change of intensity of the surface 
low center. In addition, in limited areas where the direc- 
tion of motion of the radiation Low is available, the upper 
air pattern and 500-mb. winds may be approximated. 
A special family, is 
recognizable and the general location of the associated 
maximum winds may be determined. In all latitudes the 
radiation analysis whether determined by radiation data 


synoptic situation, the eyeclone 


alone, nephanalysis alone, or a combination of both, shows 


good continuity (figs. lla and b). 

It may also be possible to determine vertical motion 
patterns from the radiation map patterns, particularly 
when combining the low resolution, quantitative radiation 
map figures from Explorer VII with the high resolution, 
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qualitative cloud pictures from TIROS. Research on 

this relationship is being conducted at the University of 

Wisconsin. 

Radiation measurements from two satellites in polar 
orbits would increase observation density and make the 
maps more nearly synoptic, increasing their accuracy and 
value. Data storage capability aboard a single satellite 
would make it possible to have radiation maps for the 
entire earth using just single, or a few, readout stations. 

The conclusions arrived at are based on data not always 
as complete as desirable and on maps which are both non- 
synoptic and few in number. However, these limitations 
are mitigated by the consistency of the relationships over 
the 3-day period and, in the northern latitudes, by the 
density of radiation observations which lessened subjectiv- 
ity. 
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ABSTRACT 


TIROS I satellite pictures for a 3-day period about a developing ‘“‘Low system” over the central United States 


are shown. 
(classical) distribution. 


subsidence and advection of dry air along a low-level wind maximum 


Cloud distribution about surface fronts and upper-air troughs was examined and compared to the “ideal” 
The pictures and analyses reveal that cloud patterns were significantly changed by strong 


The cloud patterns when compared to the 


600-mb. vertical motions computed routinely by the Joint Numerical Weather Prediction Unit showed only partial 


agreement. 


A strong relationship between low cloud cover and surface relative humidity was found; this relationship may be 
useful under certain conditions in determining the stability of the air over areas where upper-air data are sparse. 


1. INTRODUCTION 


This study was undertaken to investigate the cloud 
distribution, as indicated by satellite pictures in addition 
to the usual synoptic methods, for several stages in a 
developing low pressure system in the central United 
States. Analvses of the pictures and meteorological data 
were made in the following manner: 

1. Latitude-longitude grids were drawn and superim- 
posed on selected satellite pictures taken over the low 
pressure system on May 19 through May 21, 1960. 

2. To the pictures were added surface analyses and 


selected data. For example see figure 1. 


3. Surface analyses and upper-air analyses were made. 
(Only the surface maps are shown with the 500-mb. 
trough positions superimposed.) For example — see 
figure 3b. 

4. Simplified schematic diagrams made from cloud 
pictures were transferred to polar stereographic maps of 
Over the United 
States cloud types were taken from the nearest 3-hourly 


1800 or 2100 GMT). 


cloud data and frontal positions. 


observation time (i.e., For example 
see figure 3a. 

5. The relationships between cloud distribution and 
(a) synoptic situation, (b) low-level wind maximum, 
¢) vertical motions, (d) surface relative humidity, were 
examined for each of the three days. 


2. CLOUD DISTRIBUTIONS AND ANALYSES, 
MAY 19, 1960 


Figures 1 and 2 picture the cloud fields which are 
schematically represented on figure 3a. Figure 3b is the 
corresponding surface analysis with the 500-mb. trough 


position superimposed. 


research has been supported by the Nationa Aeronautics and Space Adminis 


SYNOPTIC SITUATION 


Figure 3a shows that the bulk of the cloudiness associ- 
ated with the surface cold front lay well to the rear of the 
front and to the northwest of the wave centered in 
southern Kansas. Further, the region immediately about 
and ahead of the cold front was clear, with the exception 
of the small zones of scattered to broken cumulus clouds. 
This type of distribution differs completely from the 
“ideal” distribution model [1] in which the bulk of the 
cloudiness lies ahead of the cold front with rapid clearing 
to the rear. 

The bulk of the cloudiness to the rear of the cold front 
was located in and ahead of the 500-mb. trough and agreed 
with the findings of Oliver and Oliver [2] that clouds 
are located at and ahead of their respective upper-air 
troughs. There were small zones of cloudiness behind the 
500-mb. trough; these however, were probably the result 
of orography since they lay over the highest peaks of 
the Rockies. 

The brightness and high reflectivity of the homoge- 
neous-appearing clouds immediately ahead of the warm 
front (fig. 2) can probably be attributed to the high water 
content of the clouds since continuous precipitation was 
This low cloud mass, whose trailing 
front, 


reported in this area. 
edge coincided closely with the surface warm 
extended about 300 miles ahead of the front where the 
cloud deck abruptly rose to 10,000 ft. or higher and 
This cloud 
distribution conformed nicely with the “ideal” warm 
frontal model [1]. 


extended about 600 mi. ahead of the front. 


LOW-LEVEL WIND MAXIMUM 
At 1800 Gut, May 19, a low-level wind maximum at 
5,000 ft. was located near Louisville, Ky. and extended 
into upper Michigan. A weaker wind maximum (about 
30-40 kt.) also at 5,000 ft. extended through central 
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PiGURE 1 


Oklahoma and = southern 


At this time however, there was no apparent 


Texas and northward = into 
Kansas 
relationship between the cloud distribution and the wind 
PaNUUIN 
VERTICAL MOTIONS 

The 600-mb. vertical motions, as computed by the 
Joint Numerical Weather Prediction Unit (JNWP) and 
interpolated linearly for 2000 GMT, May 19, showed that 





TIROS picture for 2001 Gur, May 19, 1960 with data and latitude-longitude grids superimposed 


clouds were found in areas of computed downward 10- 
tion and that clear to scattered regions were located in 
The magnitude of 
these motions was small For example, in the area of the 
“square cloud”’ (fig. 1 at 34° N., 99° W. 
to have been a tornado-producing cloud [3], the computed 
The comparisol of 


areas of computed upward motion. 
), which is believed 


upward motion was only 1 em. sec.~! 
the computed vertical motions with cloud pictures in this 
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FIGURE 2 


case nav not be valid due to the fact that vertical motions 
computed with the present JINWP model reflect only 
lara scale motion, while cloud fields are frequently pro- 
duces by mesoscale systems and convective overturning. 
SURFACE RELATIVE HUMIDITY 

D ing the course of this study, it was noticed that the 
surface relative humidity patterns were well correlated 
wit! Although the low clouds 
on al! three days largely reflected the influence of after- 


loud cover below 5.000 ft. 


TIROS picture for 2003 Gmt, May 19, 1960 with data and latitude-longitude grids superimposed 


noon convection, some low ceilings occurred with lronts 
where the same humidity relationship was evident, though 
less marked. 

A reasonable explanation may be that in areas where 
surface heating produces convection, a mixed subcloud 
laver should display the same humidity pattern at the 
surface and at cloud level Even under frontal clouds, 
mixing may occur in the subcloud layer, but it is likely 


that higher humidities might be produced by the cloud 
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FIGURE 3 a 
hatching and stippling, respectively. 
2000 GMT, May 19. 


(percent), 2000 Gur, May 19, 1960. 


cover itself by holding down temperatures during the 
afternoon, thus preventing the afternoon decrease of 
relative humidity that accompanies heating with a con- 
stant specific humidity. In addition, of course, clouds 
associated with a front may precipitate, in which case there 
should be a very good correlation between cloud and 
humidity, 

Some practical use might be made of this relationship, 
at least over convective areas, because it appears that the 
agreement or disagreement of surface humidity with low 


cloud ceilings depends largely on vertical mixing. This 


Schematic nephanalysis 2001-2003 amr, May 19, 1960. 
Clear, seattered, and hazy-appearing areas are left unshaded. 
(c) Broken (hatched) and overcast (stippled) cloud ceilings below 5,000 ft. with surface relative humidity isopleths 


Opaque-appearing broken and overcast regions are represented by 
(b) Surface data and analysis 


mixing, in turn, is a function of the stability; ie., where 
stability is strong surface heating can produce only few 
clouds, and where the stability is weak the heating 
quickly produces a cloud ceiling. Therefore, where high 


surface humidities can be inferred independently, lack of 


a cloud ceiling implies great stability, and a low cloud 
Even qualitative informa- 


ceiling indicates instability. 
tion of this type over data-free regions would be valuable 

Figure 3c shows the surface humidity patterns super 
imposed on the low cloud nephanalysis for 2000 Gt 


May 19. On this and subsequent figures, isopleths 0 
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FIGURE 4 


relative humidity were drawn to station data and therefore 
do not display the detail of the low cloud nephanalyses 
which were made from station reports and cloud pictures. 


\llowing for these differences, it is seen that in almost 


all of the areas where low cloud ceilings existed, the surface 


idity was greater than 70 percent. Only Belleville, 
Ill., Dayton, Ohio, McAlester, Okla., and St. Louis, Mo. 
Where the humidity 


was below 50 percent, almost no low ceilings were reported. 


¢ 


failed to report low cloud ceilings 
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TIROS picture for 1908 amt, May 20, 1960 with data and latitude-longitude grids superimposed. 


Only Waco, Tex., Cheyenne, Sheridan, and Yellowstone 
Park, Wyo. reported low ceilings in this case. 
Examination of the vertical temperature distribution 
between 1 and 2 km. at 0000 Gut, May 20, for the excep- 
stable 


over the four stations where the humidity 


tions showed that the lapse rates were quite small 
conditions 
was greater than 70 percent and no low cloud ceilings 
existed. Conversely, the lapse rates over the four sta- 
tions where low ceilings existed over low surface humidity 
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were quite large (unstable conditions) 


a curious relation- 


ship noted on the other days us well, but a relationship 
that is not of direct utility in the context just mentioned 
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> AND ANALYSES 

4 

Figures $f and 5 show the cloud patterns about the Low 
and frontal systems at 1908 and 1910 Gut, while figures 6 
and 7 show the area to the north at 2053 and 2054 G1 





The cloud schematic diagrams of these pictures were com 
bined and are shown in figure 8a. Figure 8b gives the 
corresponding surface analyses with the 500-mb. trough 
superimposed. 
SYNOPTIC SITUATI 

By 1900 amt, May 20, the surface Low in Manitoba, 
Canada had filled and a new low center had develope 
in Iowa along the cold frontal trough. This development 


marked the first strong cyclonic circulation of the | 
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Figure 6.—TIROS picture for 2053 amt, May 20, 1960 with data and latitude-longitude grids superimposed 


m. Both the warm and cold frontal systems in 
remained stationary, while the southern portion of 


cold front and the eastern portion of the warm front 


d slowly eastward and northeastward, respectively. 
very active squall line had developed and extended 
e way from Indianapolis, Ind. to Galveston, Tex. 
ow-level southwesterly winds which extended north- 
through central Texas during the development of 
juall line on the night of Mav 19 (between 0000 and 


0300 GMT, May 20) and the presence of a zone of low-level 
stability (small lapse rates) offer strong support to the 
theories of Rossby [4] and Cahn [5] and the work of 
Tepper |6] in which he speculates that ““The development 
of eastward propagating pressure-jump lines (associated 
with squall lines) was due to the addition of momentum 
to the low level southerly flow by a marked increase in the 
southerly component of winds during the night.” 

The cloud field had become organized into a broad 
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FIGURE 7 


north-south band (ig. 8a). The cloud distribution which 
lay well behind the surface cold front on the previous day 
had now assumed a position over the front and the squall 
line ina pattern following the “‘ideal’’ cold frontal model. 
Middle and high clouds still remained at and ahead of the 
500-mb. trough; but by the steepening of the surface to 
500-mb 180 to 1/50, middle and 
higher clouds were brought over the low cloud deck thus 


trough slope from 1 


forming the broad north-south band with rapid clearing 
to the rear 
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TIROS picture for 2054 Gaur, May 20, 1960 with data and latitude-longitude grids superimposed 


The breaks in the overcast band (fig. 8a) indicate the 
separation between that portion of the cloud band asso- 
ciated with the cold front and that associated with the 
squall line. These small clearing areas between 92° and 
93° W. and the weakening of the cloud band at about 
35° N. are best viewed on figure 4. It was in this region, 
at this time, that the clear tongue photographed on the 
next day (May 21) was beginning to form. 

On figures 5 and 7 note that the cloud field ahead of the 
warm front (now stationary) had become less extensive, 
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Ficgi re 8.—(a) Schematic nephanalysis 1908-1910 and 2053-2054 amt, May 20, 1960. Broken line shows axis of maximum wind at 5,000 
(from map ec). (b) Surface data and analyses 1900 and 2100 Gmut, May 20, 1960. (c) Wind data, axis of maximum wind, and_iso- 
ichs (kt.) at 5,000 ft., 1800 amr, May 20, 1960. (d) Broken and overcast cloud ceilings below 5,000 ft. with surface relative humid- 

ity isopleths (percent) 1900 and 2100 Gut, May 20, 1960. 
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850-mb. horizontal moisture advection (gm. kg.~! hr.~!), 


0000 cmt, May 21, 1960 


FIGURE 














Lapse rates between 1 and 2 km. (°C. km.~!), 0000 


GMT, May 21, 1960 


Figure 10 


less dense, and somewhat disorganized from that seen on 
the previous day. The zone of scattered clouds extending 
across and ahead of the front in the vicinity of 85° W. is 
believed to have been caused by horizontal advection of 


dry air and is discussed below. 
LOW-LEVEL WIND MAXIMUM 


Figure 8c shows the wind analysis at 1800 Gat, May 20, 
1960 for the 5,000-ft. level; the axis is also shown on 


figure 8a. The clearing in the region of the maximum wind 
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and its long narrow shape suggest it may have been pro- 
duced by the advection of dry air in the low-level jet. 

To examine this possibility, the horizontal advection of 
moisture was computed for 0000 Gmt, May 21, and is 
shown in figure 9. This computation revealed that the 
strongest advection was taking place south of the wind 
maximum. The advection of dry air across the stationary 
front near 85° W. was primarily the result of a large 
moisture gradient while the advection to the south of the 
jet maximum was caused by a large moisture gradient as 
well as high wind speeds. Because these large moisture 
gradients existed within the same air mass, it appears that 
vertical motions were creating extensive pockets of dry 
air, a matter that is investigated next. 


VERTICAL MOTIONS 


The 600-mb. JNWP vertical motions, interpolated for 
picture time, showed upward motions between 0 and | 
cm. sec.~? over the eastern half of the United States with 
a very small zone of upward motion greater than 1 em. 
sec.-' centered near Jackson, Miss. The zero isopleth 
closely delineated the trailing edge of the broad north- 
south cloud band to the rear of the cold front. The 
clearing associated with the low-level wind maximum, the 
zone of scattered clouds across the stationary front in the 
vicinity of 85° W., and the region of clear to scattered 
clouds ahead of the s juall line (see fig. 8a) were all located 
in the region of computed large-scale upward motion, 


SURFACE RELATIVE HUMIDITY 

Figure 8d shows the 1900 and 2100 Gat surface humidity 
patterns superimposed on the corresponding low cloud 
nephanalysis for May 20, 1960. The 70 percent humidity 
isopleth again showed a tendency to outline the regions of 
low cloud ceilings while there was generally an absence of 
low cloud cover where the humidity was below 50 percent. 

Of the 215 stations shown, 41 violated the 
stipulated 50 percent and 70 percent limitations: 8 
stations, located in areas labeled “DRY” showed broken 
low cloud ceilings with surface humidities less than 50 


stations 


percent; 33 stations, located in areas labeled ‘‘ MOIST”, 
had humidities greater than 70 percent but no low cloud 
ceilings. 

Examination of the lapse rates about the 850-mb. level 
approximately 5 hours after picture time indicated that 
the violations were again related to the low-level stability. 
In figure 10, which shows the lapse rates between 1 and 
2 km. for 0000 amt, May 21, it is seen that the lapse 
rates over the areas corresponding to the regions labeled 
“MOIST” were quite stable; while over the regions 
labeled “DRY” the lapse rates were quite unstable. 


4. CLOUD DISTRIBUTIONS AND ANALYSES, 
MAY 21, 1960 


Figures 11 and 12 show the spiral cloud band about 
the cold front and the clear tongue to the rear at 1816 
and 1817 cmt. Figures 13 and 14 show the northern por- 
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Figure 11.—TIROS picture for 1816 Gut, May 21, 1960 with data and latitude-longitude grids superimposed 


ion of the clear tongue and the spiral band associated with 
vortex at 2000 and 2001 emt. The combined sche- 
tic nephanalysis for these photos is shown in figure 15a. 
ure 15b gives the corresponding surface analyses. The 
-mb. trough is not shown because there was uniform 
vature about the closed system with no pronounced 


rough. 


SYNOPTIC SITUATION 


By 1800 Gmnrt, May 21, the surface low center had 


moved 180 mi. northeastward while the upper-air Low 


had deepened, become clesed, and moved approximately 
300 mi. to the northeast; this movement made the axis 
of the low system very nearly vertical. During the 
previous 23 hours, the winds in the low levels had under- 
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Figure 12.--TIROS picture for 1817 Gut, May 21, 1960 with data and latitude-longitude grids superimposed 





gone a 10-kt. increase and the rapid movement of the A large change in the cloud distribution had taken 
cold front produced a long bent-back occlusion over the place from that shown in figure 8a for the previous day. 
Great Lakes. The cold fron‘, especially south of 35° N., The cloud field associated with the cold front had become 
was undergoing frontolysis and had become shallow and a very narrow band lying over the upslope side of the 


diffuse. Appalachians and an expansive clear area with a long 
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Fiagure 13 


narrow curved clear tongue had developed behind the 
front. The broken section of clouds about 200 mi. east 
of the cold front partly reflected the weak instability line. 

Cloud patterns on this day complied very nicely with 
the “‘ideal”’ distribution model. 
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TIROS picture for 2000 Gur, May 21, 1960 with data and latitude-longitude grids superimposed 


W-LEVEL WIND MAXIMUM 


Figure 15c shows the wind analysis at 1800 Gut, May 


21, 1960 for the 5,000-ft. level. Figure 15a shows the 


nephanalysis with the axis superimposed. Notice that 
the axis coincided very closely with the sharp edge of the 
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Figure 14.—TIROS picture for 2001 gmt, May 21, 1960 with data and latitude-longitude grids superimposed 


low cloud mass which was associated with the vortex to the 
north. 

Examination of the vertical distribution of winds in this 
area showed that the jet stream at 35,000 ft. was posi- 
tioned very nearly vertically above the maximum at 5,000 
ft., while the intermediate axes about the 14,000-ft. level 
were located approximately 100 mi. toward the south. 

The horizontal moisture advection computed for 1200 
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GMT, May 21 is shown in figure 16. Again, the strongest 


advection of dry air was occurring to the south of the 


maximum axis. 


VERTICAL MOTIONS 


The 600-mb. JNWP vertical motions showed broadseale 


downward motion over the Plains States west of the Mis- 


sissippi Valley and a broad zone of upward motions 














MONTHLY WEATHER REVIEW 


NoVEMBER 1961 





< 
x 


> WZ. 
ae 


+ 
_0ererens 
econ 





























2 T ° i * 
‘ i 7 * 

\ w, d : ; . 
\3 2 | _._—_+- fo ' a 
Broken line shows axis of maximum wind at 5,000 
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relative humidity 











(a) Schematic nephanalysis 1816-1817 and 2000-2001 Gut, May 21, 1960. 


Figure 15. 


ft. (from ¢). 
kt.) at 5,000 ft., 1800 amt, May 21, 1960. 


isopleths (percent), 1800 and 2000 amt, May 21, 1960. 


(b) Surface data and analyses 1800 and 2000 gmt, May 21, 1960. 
(d) Broken and overcast cloud ceilings below 5,000 ft. with surface 
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horizontal moisture advection (gm. kg.~! 


1200 cmt, May 21, 1960. 


850-mb. 
hr.~'), 


Figure 16 


gulfing the clear tongue and a portion of the large clear 


area to the rear of the cold front. 

Careful nephanalyses of several sets of surface observa- 
tions indicated that the strongest development of the clear 
tongue occurred between 1800 cmt, May 20 and 1200 Gur, 
May 21. For that reason, the mean vertical motions for 
the period 0000 to 1200 Gur, May 21 were computed by 
the adiabatic method [7] for the 850-mb. and 700-mb. 
levels. Figure 17 shows the vertical motions at 850 mb. 

Areas near the Great Lakes were covered with clouds so 
that the computed vertical motions were based on moist 
adiabatic processes. At the grid points on figure 17 show- 
ing no values, the actual lapse rates were approximately 
equal to the moist adiabatic lapse rates and hence these 
computations failed because the stability approached zero. 
While it is not entirely satisfactory to use the adiabatic 
method at the 850-mb. level because of the influence of 
surface radiation, it is very likely that the downward mo- 
tion vielded by this computation was correct in sign. This 
follows because the diabatie cocling of air from 7:00 p.m. 
to 7:00 a.m. would produce fictitious upward motions in the 
computation, so that downward motions were shown de- 
Upward motions, by 
the same token, are suspect and the vertical motion over 


spite the masking effect of cooling. 


the Appalachian Mountains may be spurious because the 
diabatic influence was maximized in that region. 
Reliability of the computed motions was further at- 
tested to by the lapse rates about the 850-mb. level. Com- 
parison of vertical motions at 850 mb. (fig. 17) and the 
lapse rates between 1 and 2 km. (fig. 10) indicated that the 


tongue of small lapse rates, about 3°-4° C. km.~', was 


Figure 17.—850-mb. vertical motions (em. sec. ~!) computed by the 
adiabatic method, 0000-1200 amr, May 21, 1960. 


being produced by the similarly oriented tongue of down- 
ward motions. 

Vertical motions at the 700-mb. level were similar to 
those computed at 850 mb., with some slope toward the 
north of the center at 35° N., 95° W. The secondary cen- 
ter of downward motion at 40° N., 90° W. was overlain 
with upward motion at 700 mb. In summary, from these 
computations it appears that subsidence in the low trop- 
osphere was concentrated on the anticyclonic shear side 
of the low-level jet and was producing dry air which was 
then advected rapidly into the long clear tongue that is 
shown in the pictures on May 21, 1960. 


SURFACE RELATIVE HUMIDITY 


Figure 15d shows the surface relative humidity patterns 
and the low cloud nephanalysis for May 21, 1960. The 
70 percent humidity isopleth again closely delineated the 
areas of low cloud ceilings. Of the 235 stations for which 
surface humidity was available, not one station reported 
a low cloud ceiling when the humidity was below 50 per- 
cent; only 5 stations violated the prescribed 70 percent 
limitation by showing no low cloud ceilings with greater 
than 70 percent relative humidity. 

Over these 5 stations, the lapse rates at 0000 Gw7, 
May 22 again indicated that low-level stability was pre- 
vailing in the regions of the violations. The increase in 
area of greater lapse rates (greater than 6°C.km.~') from 
that of the previous day further suggests that surface 
humidity is better correlated with low cloud cover when 
the air is nearly unstable and substantial mixing is occur- 


ring. 
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5. SUMMARY 


Over the 3-day period, cloud distribution about the sur- 
face fronts conformed nicely with the “ideal” distribution 
models with the exception of the distribution about the cold 
front on May 19, 1960. On this day, clouds were found 
well to the rear of the cold front because the 500-mb. 
trough, with which the bulk of the clouds was associated, 
was located approximately 500 mi. to the rear of the sur- 
face frontal position. In such north-south cold frontal 
systems, the amount of cloudiness behind the surface 
front may be an indicator of the slope of the upper-air 
system. 

The long clear tongue shown in figures 12 and 14 
developed in the region of a low-level wind maximum. 
This striking feature, which appears to be a frequent 
mark of this type of synoptic situation, was caused by the 
development of a strong band of subsidence on the anti- 
cyclonic shear side of the low-level wind maximum, This 
strong subsidence at both the 850- and 700-mb. levels, on 
the anticyclonic shear side of the maximum wind speed 
intensified the horizontal advection of dry air and hence 
the clearing far downstream from the isotach maximum. 

The 600-mb. JNWP vertical motions computed by the 
present model reflected some large-scale vertical motion 
fields. 
computed downward motion and often no clouds were 
The vertical motions 


Clouds were sometimes photographed in areas of 


found in areas of upward motion. 
computed by the adiabatic method for the period 0000 
1200 amt, May 21 agreed quite well with the cloudiness 
observed at 0600 cmt, May 21. 

There was a high correlation between surface humidity 
and low cloud cover. Over the 3-day period, of the 660 
stations examined, only 54 violated the 50 percent and 
70 percent limitations and each of these 54 violations fell 
within 120 miles of the prescribed humidity isopleth. In 
general, the following relationships were seen: (a) low 
cloud cover was Overcast where surface humidity was 80 
percent or greater; (b) low cloud ceilings (broken to over- 
cast) were observed where the surface relative humidity 
was 70 percent or greater; (c) relative humidities of 50 
percent to 70 percent were the middle ground where low 
cloud coverage was both smaller and greater than 5/10; 
(d) low cloud ceilings were not reported where surface hu- 
midity was below 50 percent. 
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The strong relationship between suiface humidity and 
low cloud cover appears to be characteristic of summer 
afternoon regimes where strong surface heating produces 
instability in the low levels and thus allows substantial 
mixing of the subcloud layer. It was seen that the poorest 
agreement of the surface humidity and low cloud cover 
existed over areas of very small lapse rates (stable con- 
ditions). On the other hand, when lapse rates were larger 
(more unstable conditions) there was considerable ver- 
tical mixing of a moist layer such that surface humidities 
were closely related to low cloud cover. 
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ABSTRACT 


On May 16, 1960, TIROS [ televised a series of pictures over the central United States during a severe local storm 


situation. 


Comparisons are made between the location of the isolated, bright cloud elements seen in those pictures 


and (1) the location of surface weather and synoptic features, (2) the location of sferics observations, and (3) the 


location of the upper-level jet maxima believed important in the development of severe local storm activity. 
parison is made of the cloud cover shown in the pictures to the cover indicated by the surface observations. 


Com- 
Apparent 


brightness and cloud cover of the pictures is discussed in relation to both the weather activity occurring at the surface, 


and cloud types and cover observed at the surface. 


Also, cloud streets which appear in the photographs are dis- 


cussed relative to the observed wind patterns and vertical temperature structure. 


Evidence is presented indicating that the severe local storm activity occurred beneath three bright mesoscale 


cloud masses 


scale patterns. 


the more vigorous activity occurring beneath two masses which were largely isolated from the large- 
Further, the latter two clouds masses are shown to be correlated well in location with sferics observa- 


tions and jet stream maxima in the lower levels of the atmosphere 


1. INTRODUCTION 


bright cloud 
appeared in the many pictures received from the TIROS 
limited of 
pictures have shown small-scale very bright cloud patches 
These 


patches, because of their relative isolation, strong reflec- 


Large-scale, areas have frequently 


meteorological satellites. A number these 


somewhat isolated from the large-scale patterns. 


tivity, and mesoscale proportions, have suggested intense 
convective activity, such as a system of well developed 
thunderstorms perhaps capable of producing a family of 
In at 
such a case has been documented [11]. 

On May 16, 1960, during pass 659, TIROS I televised 
a sequence of pictures showing bright, isolated masses of 


severe local storms. least one instance to date, 


clouds over the central United States during a severe 
weather situation. A study of that case was undertaken 
to provide evidence supporting the hypothesis that these 
distinctive cloud masses were associated with the occur- 


rence of severe weather in the Middle West. In particular, 


these clouds were investigated relative to synoptic ob- 


In addition, the 
photographs were compared with surface cloud observa- 
Also, 
rather prominent wave clouds or cloud streets appearing 
in the satellite pictures were examined relative to the 


servations, analyses, and _ sferics. 


tions to aid in identifying the clouds in the pictures. 


theoretical models of cloud streets. 


is research has been supported by the National Aeronautics and Space 


nistration 


2. DISCUSSION OF PICTURES 


Five wide-angle TIROS I pictures were photographed 
at approximately 30-second intervals between 205433 
GMT and 205831 cmt on May 16, 1960 over the central 
United States. 
State borders have been superimposed on each of the five 


To facilitate discussion of the pictures, 
pictures (figs. 1-5). Rectification and gridding was done 
using the procedures discussed in [10]. 

Two of the three cloud features of interest (A, B, and 
C) are shown in the figures 2-4. Note that the three 
major bright areas were located (A) in northern Missouri 
and southeastern Iowa, (B) over the Illinois and Missouri 
border, and (C) over Indiana. The cloud mass C, since 
it is not isolated as are the clouds A and B, does not stand 
out as clearly. While the western border is well defined, 
very bright, and easily distinguished, the eastern border 
blends into the large-scale cover. 

The cloud masses described above appeared consider- 
ably brighter on the transparencies than in reproduction 
here. Unfortunately, prints of satellite pictures do not 
usually differentiate scales of brightness as adequately as 
do the original 35-mm. transparencies. 

Several other prominent and interesting cloud patterns 
appearing in the pictures are worthy of mention here. 
For instance, over central Oklahoma and northern Texas 
(figs. 1 and 2) the clouds are organized in parallel lines 
or streets. The lines are rather closely spaced and appear 
comprised of small cloud elements much like cumulus 





Figure | 
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Gridded photograph of frame 4 from the TIROS I pass 659, May 16, 1960, 2100 Gm1 


Cloud streets are seen at the right 


in the picture 


Farther to the east, a long, moderately brigl.t 


cells. 
cloud band may be seen oriented northeast-southwest 
through central Arkansas (fig. 2). This band suggests a 
line of cumuliform clouds. 

In figure 3 the moderately bright, extensive cloud cover 
across the top of the picture was located in the western 
Lakes area westward into Minnesota and northwestern 
Iowa. Immediately to the east, over the eastern Great 
Lakes area (fig. 5), the clouds became filmy and streaked, 
suggesting that they were cirriform. 


Figure 6 is a schematic of the essential features in the 


five TIROS I features are 
described by the prominence of the boundaries and relative 
brightness. 


photographs. The cloud 


This presentation, hopefully, will permit the 
reader to view the features in composite. 


3. COMPARISONS OF PICTURES AND METEORO:- 
LOGICAL ANALYSES 


SURFACE ANALYSIS 


The surface chart prepared from the 2100 Gur surface 
airways Observations corresponds to the picture times. 
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Gridded photograph of frame 7 from the TIROS I pass 659, May 16, 1960, 2100 Gm1 


Two of three very bright cloud masses 


associated with severe weather are designated A and B 


A series of hourly charts was prepared to observe short- 
period changes and to check the continuity of the 2100 
GMT analysis. In order that these analyses might be as 
objective as possible, they were performed independently 
of the picture rectifications (geographical location of the 
cloud elements). 

In the interest of brevity, 


y, only the 2100 Gar surface 
7). Two small waves appeared 


analysis is shown here (fig. 
on the cold front—one in eastern Kansas and another in 
western Oklahoma. <A front from the 


warm extended 


pressure center southeastward into West 
Virginia. Continuous light rain fell to the north of the 
lowa low center. A few were 
reported in Michigan north of the warm front. 

Of primary interest were the thunderstorms and associ- 
One 
group of thunderstorms produced some hail along the 


group 


lowa low 


pressure rain showers 


ated severe weather reported in the warm sector 


cold front in northern Missouri; a second was 


accompanied by tornadoes and strong winds near the 
Illinois-Missouri border; and a third group spawned a 





450 





> 


FIGURE 3 


same two which appear in figure 2. These two masses 


lone tornado in Indiana. The convective systems are 
designated A, B, and C, respectively, in figure 7. 

Severe weather occurred over a span of 5 hours between 
1830 and 2330 omr (fig. 8) but a good portion of the 
activity was observed within one hour of 2100 Gmr—the 
time of the TIROS pictures. Though the severe weather 
during this period was outside the Severe Local Storms 
Research Network, four mesoscale pressure systems could 
be detected (fig. 7)—two mesohighs and two mesolows. 


These systems were analy zed from the airways observa- 
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Gridded photograph of frame 11 from the TIROS I pass 659, May 16, 


are 
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1960, 2100 G1 
isolated from the 


A and B are the 


large-scale cover to the north. 


The cloud masses 
somewhat 


tions through close examination of the behavior of the 
The 


analysis of these pressure systems is crude in comparison 


wind, pressure, temperature, and weather in time. 


with that which might be prepared from a dense network 
of accelerated microbarograph and thermograph traces 
[5).2 
sufficient evidence that pressure disturbances of the type 
indicated did exist in 


For the scope of this report there appears to be 


the areas shown. In short, the 


2 Dr, Tetsuya Fujita of the University of Chicago has been performing a very detuiled 
mesoanalysis of this situation 
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Gridded photograph of frame 15 from the TIROS I pass 659, May 16, 1960, 2100 Gur. The third bright cloud mass appear 
here as area C. This cloud mass is not isolated from the large-scale cloud cover as are clouds A and B 


severe weather activity in eastern Missouri and that in The only tornado associated with it developed 45 minutes 
Indiana were each associated with a mesohigh and the earlier. Its history dated back several hours, having 
accompanying pressure rise or pressure jump line. It earlier passed through central Illinois and northern 
appears that the storms in northern Missouri were Missouri. In all likelihood the system was in a state of 
associated with the cold front, although one should not decay. 

rule out the possibility of a mesosystem there. The There was no definite indication of unusual weather 
location of the cold front in Missouri was somewhat activity within the two mesolows. Possibly, though, the 


nebulous because of the convective activity present. mesolow in northern Indiana was a tornado-Low or the 


lhe convective system C, over Indiana, was less severe remnants of the tornado observed in Indiana (figs. 7 


at 2100 Gut than were the two disturbances to the west. and 8). 
613342—61——5 
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hIiGURB 5 
cover appears at the left 
the Great Lakes area. 


The weather within 1 hour of 2100 
amr, and the major surface synoptic features and the 
the 2100 


shown in composite with the rectified cloud schematic 


severe reports 


surface weather from GMT surface chart are 


in figure 9 to expedite comparisons. One notes immedi- 
ately the excellent the location of 
the thunderstorm and the severe storm reports relative 


correlation between 


to two of the three bright cloud masses—one in northern 
Missouri and southeastern Iowa (area A), and the other 
in eastern Missouri and southwestern Illinois (area B). 


Gridded photograph of frame 19 from the TIROS I pass 659, May 16, 
The filmy and streaked clouds through the center of the picture are associated with cirrus reported In 
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ar ae 
\ 


1960, 2100 amr. A portion of the large-scale cloud 


The cold front “slices’’ through the center of the first mass 
and a squall or pressure jump line through the second. 
Geographical agreement between the mesosystem and 
the third bright cloud mass, C, 
the two the west. Thunderstorms 
reported both within and without the bright cloud while 
the squall line, as analyzed, is not aligned with the long 
dimension of the bright cloud. Lack of agreement may 
have been influenced by two factors. One, the eastern 
portion of cloud mass C was very diffuse with the bright- 


is not as striking as for 


systems to were 
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Figure 6.—Schematic rectification of the cloud features appearing in figures 1-5. 
storm systems in figure 7. 
Superimposed on the schematic are the 2100 emt sferics locations taken from figure 10. 


cally lettered cloud masses in figures 2-4 and the severe 
subjectively from the pictures. 


ness decreasing toward the east, eventually blending 
into the large-scale cloud cover (fig. 4) and thus precluding 
Two, since the 
disturbance in Indiana was old and decaying, it was not 
sharply defined. 


clear delineation of the eastern border. 


SFERICS 

Figure 10 represents the 50 ke. sferics observations 
from four station fixes during the interval 2055-2100 
eur, May 16, 1960. Quite obviously there was distinct 
and frequent sferics return from the area of northern 
and northwestern Missouri and southeastern Iowa and 
from the area of east-central Missouri and southwestern 
Ilinois. 
sferics reports correspond extremely well in location with 
two of the three areas of small-scale convective activity, 
areas A and B of figure 6, described earlier. In fact, these 
sferics reports probably outline, in greater detail than the 


Comparison of figures 7 and 10 shows these 


Areas A, B, and C correspond both to those identi- 
All brightness descriptions were determined 


surface observations, those areas where rather vigorous 
convective activity was taking place. 

There was little or no 50 ke. sferies return of conse- 
quence at this time in Indiana, area C, even though 
observations of thunderstorms were made there. Several 
explanations may be given for this but a plausible one is 
that the Indiana storms, being less active and decaying, 
were emitting electrical disturbances in a different 
frequency range [4] than those to the west and southwest 
where the systems were vigorous and developing or 
regenerating, producing severe storms until 2330 emr. 

Since two convective areas (A and B) were shown to 
have been within the bounds of bright clouds seen from 
TIROS I, then the sferics observations and bright clouds 
should be expected to exhibit a measure of correlation. 
Figure 6, a composite of the sferics return and cloud 
schematic, reveals a correlation beyond expectation. The 
bright cloud areas match the greater portion of the sferics 
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FIGURE 7. 
observed cloud types, and observed weather. 
sive categories are lighter in shading. Category 1 


scattered) is unshaded Cloud types, 


applies 


return area to the extent that the boundaries are coincident 


at many points, 
UPPER AIR STREAMLINE ANALYSES 


Streamline analyses were performed (after the technique 
of Palmer et al. [9]) for the 5,000- and 20,000-ft. levels 
from the 1800 Gmr May 16, 1960 and the 0000 Gur 
May 17, 1960 upper air observations over the central 
portion of the United States (figs. 11 12). The 
object of the analyses was to determine the configuration 


of low- and high-level jet streams and jet maxima for 
the purpose of locating a convective mechanism, as 


proposed by Beebe and Bates [1], which might be related 
to the bright clouds appearing in the TIROS I pictures 


from pass 659. 

In short, the Beebe and Bates divergence models for 
location of the convective mechanism were evolved from 
a determination of the divergence fields about jet maxima 


The cloud cover is separated into the four categories described in the text. 
essentially overcast 
weather, fronts, ete., are 


NOVEMBER 


SURFACE ANALYSIS 
2100 GMT 
I6 MAY 1960 


Surface chart for 2100 Gur, May 16, 1960, showing, in addition to the pressure and frontal analysis, the observed cloud cover, 


succes- 
is described by heaviest shading, while category 4 (clear or 


represented by the usual symbols. The legend of figure 9 also 


at both low and high levels. The divergence field may 
be determined from the divergence term of the vorticity 
equation by evaluating from the wind field alone the 
advection of relative vorticity about a jet maximum. 
If the configuration of a high-level jet maximum over a 
low-level jet maximum is such that divergence is occurring 
over convergence, Beebe and Bates propose that the ver- 
tical motion and modification of the temperature distri- 
bution would release convective instability. 

The analyses shown in figures 11 and 12 cannot be 
considered unique because of the sparseness of upper-air 
observations. A quantitative evaluation of the diver- 
gence term from these analyses would then be somewhat 
presumptuous for the scale we are considering vere. 
Therefore, the models were applied on the 5,000- and 
20,000-ft. the 
at the times considered. The models are the cases of 


charts to locate convective mechanisms 


(1) a jet stream of no curvature (fig. 11a), (2) a eyelonically 
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Figure 8.—Severe storm reports in the Middle West on May 16, 






curved jet with the wind maximum located at the point 
of maximum curvature (figs. 12 a and b), and (3) an 
anticyclonically curved jet with the wind maximum 
located at the point of maximum curvature (fig. 11b). 
These models are simplified compared with the many jet 
structures which may occur, but nonetheless they are 
useful tools in synoptic meteorology. 

Based on the 5,000-ft. and 20,000-ft. streamline charts 
for 1800 Gar, it was possible using the models to delineate 
with confidence one area where divergence occurred over 
convergence. That area (shaded in figs. 11 a and b) 
east of the jet axes lay over southern Missouri and north- 
ern Arkansas or just southeast of cloud mass B where 
severe weather began breaking out at 2000 Gar. 

By 0000 amr of May 17, 1960, the jet convective 


portions of southeastern Iowa, 


mechanism lay over 
northeastern Missouri, and western Indiana (shaded in 
figs. 12 a and b). Wind damage and tornadoes occurred 
in that area with system A during the 2-hour period prior 


to 1000 Gmr. Numerous thunderstorms and some hail 





1960, Greenwich time. 
homa, but these reports are not shown since they had no direct relation to the storms in the Middle West. 
















and Okla- 





Other storms occurred in Kansas 


within cloud mass A (fig. 3) were forerunners of the 
tornadoes. 

Thus, over the 6-hour interval between 1800 and 0000 
GMT, low- and high-level jet maxima and consequently 
the dynamic convective mechanism were restricted to 
the immediate vicinity of the two areas of severe weather 
associated with the cloud masses A and B. Certainly 
agreement in location between the cloud masses at 2100 
GMT and the mechanisms at 1800 and 0000 Gmr is not as 
good as it might be. One must remember that a precise 
fix of the mechanism is hampered by the wide spacing of 
the upper air observations, and, in this case, by the lack 
of time agreement of the upper air observations with the 
TIROS pictures. 

Significantly, the third group of thunderstorms (area C) 
was neither as severe as the more western systems nor 
directly under the influence of the jet maxima. During 
the interval 1800-0000 Gur, system C continually moved 
eastward away from the dynamic mechanism, lending 
further support to a decaying thunderstorm situation over 
Indiana. 
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FIGURE 9. 
and cloud types. 


SURFACE NEPHANALYSIS 

The sky condition 
according to the following four categories: 

1. Overcast low and/or middle clouds; or overcast high 
clouds with at least broken low and/or middle clouds. 

2. Broken low and/or middle clouds but high clouds 
not overcast if any appear. 

3. Broken to overcast high clouds with clear to scattered 
low and/or middle clouds. 

4. Clear skies or scattered low and/or middle clouds. 

This form of nephanalysis was adopted since it might 
better correspond to the apparent cloud cover and bright- 
ness seen in the pictures, particularly if the premise is 
accepted that high clouds are not always readily visible 
in earth satellite photographs. Experience has shown 
that cirriform clouds are dense, the TIROS 
cameras do not usually detect them. 

In the following few paragraphs, comparisons are made 


surface-observed was analyzed 


unless 


Composite of TIROS cloud schematic and surface observations. 
Other features of the diagram are given in the legend. 


Individual station plots give cloud cover, cloud heights, 


between the TIROS pictures (figs. 1-5) and the surface 


nephanalysis (fig. 7). The reader may find it beneficial 
to refer also to the TIROS schematic in figure 6. 

In the cold air the overcast area (fig. 7) as described by 
category 1 above, conforms closely to the standard model 
of heavy clouds and steady precipitation extending from 
the warm front counterclockwise around the low pressure 
center. This overcast region may be discerned in figure 3 
as an area of continuous, moderately bright cloud cover. 

A second area of overcast skies (category 1) was located 
in the northern portion of the warm sector. That are: 
contained the cloud masses A, B, and C (figs. 2 and 3). 

Areas of broken low and/or middle cloudiness (hatched 
areas of fig. 7) over Oklahoma and in the central Plains 
States are reflected in the TIROS pictures as small-scale 
patches of clouds or as individual cloud elements, sonie 
organized into cloud streets. 

Those areas on the surface nephanalysis (fig. 7) classi- 
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fied. under category 3, as broken or overcast high cloud 
(light dotted areas), are seen as either clouds of dull 
filmy appearance or areas of no cloudiness. In the 
eastern Lakes area (fig. 5), the clouds were filmy. Over 
south-central Missouri, figure 3 pictures a very dark area. 
Yet broken to overcast cirriform cloud cover with few or 
no lower clouds was reported in both locations (fig. 7). 
The observations indicated that the high clouds were 
dense over the eastern Lakes but thin over south-central 
Missouri. 


4. CLOUD INTERPRETATION 


In fic; re 9, the schematic of the five TIROS pictures 
is shown in composite with the surface cloud observations. 
From the preceding discussion of the nephanalysis in 
figure 7 and by inspection of figure 9, one may see that the 
surface observed cloud cover agrees well with the TIROS 
observed cloud cover. 

There are several positions where the rectified clouds 
are not consistent with These 
discrepancies indicate a maximum rectification error of 
about 20 miles. examples are Rantoul, Il. 
(RAN) and Terre Haute, Ind. (HUF). Both stations 
reported considerable cloudiness yet lie in a dark portion 
of the picture which indicates little or no cloudiness. 
One can see from figure 9 that only a slight geographical 
readjustment in the rectified cloud schematic would be 
required to have these stations fall within a cloudy area. 

Several interesting comparisons may be drawn from 
figure 9 with regard to observed clouds and the TIROS 
observations. For instance, stations which fell 
within the bounds of the three very bright cloud masses 
(A, B, and C) all reported an overcast cover with at least 
broken clouds in the low and middle layers. These 
stations indicated a predominance of cumuliform clouds, 
with cumulonimbus and towering cumulus or showers 


surface observations. 


Two 


those 


and thunderstorms being most prevalent. 

Observations of middle and high clouds between the 
numerous heavy cumulus structures suggest that the 
anvil tops united while spreading out under the influence 
of the strong winds in the upper levels. It is not sur- 
prising then that the cloud masses A, B, and C (figs. 3 
and 4) should exhibit continuous cloudiness rather than 
individual, large cumulus clouds. 

Figures 11b and 12b illustrate the strong flow present 
at 20,000 ft.; the flow pattern at 30,000 ft. (not shown) 
was very similar. The westerly and southwesterly 
winds over the clouds A, B, and C (fig. 9) explain the 
extension of the cloud masses to the east and northeast 
of the areas of heaviest convection. This effect is par- 
ticularly noticeable in cloud B where the cloudiness 
protrudes 50 or more miles ahead of the squall or pressure 
jump line. 

Another interesting feature of the pictures is the mod- 
erately bright cloud band oriented northeast-southwest 
through central Arkansas (fig. 2). Little Rock (LIT) 


and Texarkana (TXK) recorded towering cumulus to 
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Five-minute, four-station fixes of sferics observations 
at 2100 cmt, May 16, 1960. 


Figure 10. 


the north, and northwest through northeast, respectively 
(fig. 9). 
position some 20 miles northwest of each station. 


Both reports give credence to the cloud band 


A second apparently overcast area of moderate bright- 
ness as seen in the TIROS pictures occurred north of the 
warm front extending from the western Great Lakes are: 
counterclockwise the Iowa low pressure center 
into eastern Nebraska and western lowa (fig. 9). Station 
observations within this indicated overcast 
skies while stratocumulus was the most predominantly 


about 
area all 
reported cloud type. 


In Kansas, portions of Nebraska, Oklahoma, and north 
central Texas, the pictures (figs. 1 and 2) convey the 


impression of broken to occasionally scattered cloud 


cover comprised largely of convective cells. Many 
small, individual elements can be seen, some grouped 
together in patches of irregular size and shape, and some 
organized into cloud streets. Though the surface obser- 
vations (fig. 9) do not reveal the details of organization, 
scattered to broken cover of convective clouds was re- 
ported throughout the altocumulus castellanus 
in Texas and cumulus elsewhere. 


5. CLOUD STREETS 


area 


Cloud streets have been a source of interest 


with the advent of rocket and satellite photography. 


great 
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AO 2000 
5,000 FEET 
STREAMLINES 


1800 GMT 
16 MAY 1960 
a 


20,000 FEET 
STREAMLINES 
1800 GMT 
6 MAY 1960 
. |b 


Figure 11.—Streamline analyses for May 16, 1800 Gut: (a) 5,000- 
ft. level; (b) 20,000-ft. level. Shaded area represents the region 
under the influence of the convective mechanism as determined 
by the Beebe and Bates [1] models using the 5,000- and 20,000-ft. 
levels. Jet streams are heavy solid lines; streamlines, light solid 


lines; and isotachs, dashed lines. 


The literature has suggested that cloud streets are re- 
lated to the wind field. Meteorologists, therefore, have 
been hopeful that observations of cloud streets from 
satellite vehicles would aid in determining wind directions, 
particularly in sparse data areas. 

Observations from rockets and satellites have shown 
streets both parallel and perpendicular to the flow or 
the vertical shear. Theories supporting 
both extremes are propounded in the literature [6, 7]. 
Kuettner’s [7] theory of winds parallel to the streets 
depends on a “‘nose’’ profile of wind speed where speed 
shear changes sign with height. Haurwitz’s [6] theory, 
on the other hand, supports formation of either billow 
cloud rolls or cloud streets composed of polygonal cells 
at an internal boundary surface of density discontinuity, 


sometimes to 
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Figure 12.—Streamline analyses for May 17, 0000 Gat: (a) 5,000- 
ft. level; (b) 20,000-ft. level. See legend to figure 11. 


the rolls or streets being normal to the wind shear vector 


across the boundary. 

Surface observations and a sounding in the vicinity of 
the cloud streets observed by TIROS I (figs. 1 and 2) over 
Oklahoma on May 16 suggest that these streets conform 
to the Haurwitz theory. The cloud bases were reported 
at 4,000 to 6,000 ft. (fig. 12) while the Fort Smith, Ark. 
(FSM) sounding (fig. 13a), though taken 3 hours earlier, 
indicated a temperature inversion at approximately the 
same altitude. The wind shear vector across the inversion 
surface or internal boundary at Fort Smith was normal 
to the cloud streets as rectified. Winds immediately below 
and above the inversion were 190° at 20 kt. and 230° at 
37 kt., respectively, which determined a shear vector of 
260° at 26 kt. The cloud streets generally formed an 
angle of 80° to 90° with the shear vector. 

Haurwitz indicates that cloud rows can theoretically be 
at any angle except zero to the wind shear across an 
internal boundary surface but that normality to the wind 
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Figure 13.—(a) Fort Smith, Ark., sounding. (b) Fort Worth, 


most plausible. For the case of perpendicu- 
simplified equation for maximum wavelength 


infinite lateral extent is 


shear is the 
larity, the 
L for 


rows or streets of 


r(Al’)?T/q 
(Al’)?7 (1 —y)/gl? 


(1) 


|AT?-+ 


where g=acceleration of gravity 
AU 


[=mean absolute temperature of the two layers 


shear speed across the internal boundary 


temperature difference across the boundary 
adiabatic Japse rate or moist adiabatic rate if 
condensation occurs 
lapse rate of the atmosphere which is assumed 
equal either side of the boundary. 
When the waves are of finite lateral extent, as when the 
streets are composed of cells, the wavelength according to 


equation (1) is modified as follows: 


L 


. [1+L</a*]'/* 


where L equals wavelength of streets of finite lateral 
extent, and \ equals lateral wavelength or lateral spacing 
of cells. 

Equations (1) and (2) were applied in the May 16 cloud 
street case using values from the 1800 Gaur Fort Smith 
Though observed 3 hours earlier, the sounding 
The theoreti- 


sounding. 
Was nearest the streets in time and space. 
al wavelengths or spaces between streets thus determined 
e 3.3 miles and 1.7 miles (assuming \ to be 1 mile) for 
es of infinite and finite lateral extent, respectively. 


1. 
-3%0 -20 -10 [e) 


Tex., 
tribution; the broken line, the vertical dew-point distribution; and the light solid line, the 300° K 
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sounding. The heavy solid line is the vertical temperature dis- 


adiabat 


Since the cloud streets in the TIROS pictures appear to be 
made up of individual cloud cells, then the theoretical 
wavelength of 1.7 miles would seem to apply. However, 
the surface observations of cumulus congestus and size of 
the cloud elements in the pictures indicate the cellular 
structure is more complicated than the simple cells in 
streets of finite lateral extent. Then for purposes of this 
discussion the theoretical spacing of 3.3 miles as an initial 
condition will be considered. The cloud street spacing 
in the TIROS photographs is approximately 8 or 9 miles 

Riehl et al. [8] have noted apparent increases in cloud 
street spacing through enhancement of some streets and 
the 


Conover [3] 


suppression of others, particularly when cumulus 


became more vertically developed. has 


observed, through radar, major bands between which 
He, 


noted streets, using time lapse photography, that have 


were minor bands, more closely spaced. too, has 
disappeared or were suppressed within a pattern of streets 
Though not speaking specifically of streets, Byers and 
Braham [2] observed preferred grouping of cells during the 
Thunderstorm Project in Florida. 

Actually the theoretical treatment of Haurwitz is appli- 
cable only to clouds of small vertical development, whereas 
the streets seen in the pictures are composed of conglom- 
erates of cumulus cells of considerable vertical develop- 
ment. It is reasonable to suppose that the internal 
boundary waves provided sufficient mechanical vertical 
motion to release the latent instability for the formation 
of cumulus congestus. Further, as the cumulus developed 
some streets were damped while others were accentuated, 
The street 


thus increasing the wavelength observed 
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spacing of 8 or 9 miles and a theoretical spacing of 3.3 
miles suggest that every third street was enhanced and 
those in between suppressed. 

lhe presence of cloud streets over extreme north central 
Texas (fig. 1) is not easily explained. The clouds were 
based at 12,000 ft., much higher than the inversion at 
Fort Worth (fig. 13b), and the shear between the 12,000- 
and 14,000-ft. levels, though small, approximately paral- 
leled the streets. There is then no upparent relationship 
to the Haurwitz theory at time of the sounding. At 
cloud level, the wind and streets were perpendicular pre- 
cluding any obvious relationship to the Kuettner theory. 


6. SUMMARY 


Evidence has been presented supporting the premise 


that the isolated brightest mesoscale cloud masses seen 
in the TIROS pictures of May 16, 1960 were organized 
systems of numerous thunderstorms capable of producing 
the 


developing where the strongest convective mechanism was 


severe local storms. Further, two of masses were 
present; that is, the clouds were associated with low- and 
high-level jet stream maxima and surface fronts or pressure 
jump lines. The third bright area, though not isolated 
at picture time, was a manifestation of a decaying group 
of thunderstorms which had moved eastward for several 
hours, away from the jet convective mechanism to the 
Wwesl, 

Then, with regard to cloud interpretation, it is not 
surprising that the bright masses were concentrations of 
cumulus congestus and cumulonimbus clouds observed 
from the surface. In general, the brightest areas were 
both overcast and cumuliform in nature, the smaller-scale 


areas representing cumulus of vertical development and 
Scattered to broken 
cumulus or cumulus congestus clouds appear in the pic- 


the large-scale area, stratocumulus. 


tures as darker areas comprised of small cells or irregular 
bright patches, some organized and some not. Cirriform 
clouds, overall, gave the poorest reflectivity or none at all, 
particularly where thin. 

The cloud streets uppear related to the theory suggest- 
hetween wind shear 


ing a perpendicular orientation 


through an inversion and cloud streets which form there 
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We have seen that TIROS | has at least on two occa- 
sions, May 16 and 19, 1960, clearly demonstrated the 
capability of television-equipped meteorological satellites 
to detect severe-weather-producing cloud masses, both 
prior to storm development and during storm occurrence. 
In the future, vastly improved television coverage is 
anticipated in both time and space so that satellites are 
likely to be effective tools for both detection and tracking 
of mesoscale severe weather systems. 
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ABSTRACT 


A comparison is made between the stability properties of the non-divergent and divergent one-parameter models. 





It is shown that the introduction of divergence in the one-parameter model reduces the rate of growth of unstable 
disturbances and confines the instability to a more narrow band of wavelengths. 

Changes in zonal momentum, momentum transports, and energy conversions between mean flow kinetic energy 
and eddy kinetic energy are investigated in the linear case, as well as by extended integrations of the spectral form 


of the prognostic equations allowing only a few wave numbers. 
Long-term variations in barotropic flow, where the flow is initially stable, are investigated using the spectral 


formulation, but allowing only as many wave numbers as are needed to investigate variations in the profile of the 


zonal wind. 





1. INTRODUCTION 


The stability of barotropic flow has been investigated 
in some detail in earlier studies, notably by Kuo [8] 
und Eliasen [7]. The extensive use of the barotropic 
model in operational forecasting has caused a continuing 
interest in the performance of this model. Certain 
modifications have been introduced in the forecast equa- 
tions to alter the forecasts of planetary waves. These 
semi-empirical corrections are described by Cressman 
5} and Wiin-Nielsen [17]. The so-called divergent 
baratropic model is probably going to be used in opera- 
tional forecasting for some time to come. It is therefore 
of interest to study the stability of disturbances in this 
model and to compare the results with the corresponding 
results for the non-divergent, one-parameter model. 
The first part of this paper is concerned with such a study. 

There has been recently an increasing interest in the 
performance of the numerical prediction models when 
they are integrated over extended periods of time up to 
several days, weeks, or even months. Thompson [15] 
has developed a heuristic theory of long-period velocity 
Variations in barotropic flow. It was especially shown 
from this theory that intense jet stream maxima have 
au tendency to split into two maxima traveling toward 
the north and the south, a situation which seems to be 
characteristic of certain types of atmospheric flow con- 
nected with blocking. Lorenz [9] has shown how the 
barotropic vorticity equation can be reduced to a maxi- 
i1um simplification and has investigated the extended 


lime-integration of such a model. The same approach 


was used by Bryan [3] to perform experiments with a 
ple baroclinic model dealing with the different regimes 


Ss 
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caused by changes of the rotation of the earth and the 
rate of heating 

Saltzman [14] has investigated the non-linear inter- 
action between the zonal flow, large-scale, and small- 
scale disturbances by an extended integration of the 
barotropic vorticity equation using a limited number of 
Fourier components. 

The general circulation experiments performed by 
Phillips [12] deal with an integration of the two-parameter, 
quasi-geostrophic models over extended periods including 
adiabatic effects and friction. Extensions of Phillips’ 
work to more realistic models based on the primitive 
equations and having a greater vertical resolution in 
addition to more realistic heating and frictional effects 
are underway in different research groups. 

The second part of this paper will describe the varia- 
tions in the flow of a barotropic fluid on a large time scale. 
We shall be concerned partly with a flow which is initially 
unstable, and partly with a flow which is initially stable. 
Integral quantities and their variation in time will be given 
because we are not interested in the details of the in- 
dividual disturbances. The changes in the profile of 
the zonal wind, the momentum transports, and the 
energy conversions between the kinetic energies of the 
mean zonal flow and the disturbances will be reproduced 
in detail. We shall use a simplification of the baro- 
tropic vorticity equation for these purposes of a type 
similar to Lorenz’s minimum form, but shall allow only 
as INanvy wave components as are needed to describe a 
zonal wind profile which has at least two maxima. 


2. STABILITY STUDY OF THE DIVERGENT 
ONE-PARAMETER MODEL 






The prognostic equation for the divergent, one-param- 
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eter model will be given in the form: 


OVy/ot+V-V(VY+f)=_qdy/ot (2.1) 
where ¢ is determined by the expression : 
q _—— (2.2) 


o dpdp 


The notations have the following meaning: y is the 
stream function, VW=k Vy the horizontal non-divergent 
wind, k a vertical unit vector V?=0?/dz? +0*/d7? the two- 
dimensional Laplacian, and f the Coriolis parameter. In 
equation (2.2) o=—(a/@) 00/op is a measure of static 
stability, fp a standard value of the Coriolis parameter, and 
A (p) is the empirical function describing the vertical 
variation of the horizontal wind around the level where 
the vorticity equation is applied (500 mb.), according to 


the assumption 


Virvy. pp. O=A(pyVs(s, y, t). (2.3) 

The subscript 5 has been dropped from the start in 
equation (2.1). A derivation leading to this particular 
form of the prognostic equation has been given by Wiin- 
Nielsen [17]. 

In the first part of this section we shall consider dis- 
turbances on a zonal flow in a region bounded to the north 
and south by rigid walls. We shall divide any quantity 
into its zonal mean and deviations from the zonal means. 
For an arbitrary quantity we write: 


where 


ai adr. (2.5) 


L. is the largest possible wavelength. 
When we substitute in (2.1) and neglect second order 
terms we may derive the following equation: 


+ oF +( 8 


oT ) 4 , OY 
Oo’ or , 


or ' (2.6) 
ov" ~ Ot 

We shall consider perturbations in the eddy stream fune- 
tion of the tvpe: 


(2.7) 


yy’ (r,y,t) 


alye* 


When we substitute (2.7) into (2.6) we can write the 


latter equation in the form: 


= [ma . 1, ru, , 
(u—c) a—pairi f—;srgqe ja—0. (2.8) 
dy? dy’ 
The general eigenvalue problem to find for which 


connected values of uw and ¢ there exist solutions to 


equation (2.8), will not be treated in this paper, but we 
shall restrict ourselves in this first analytical treatment 


to an initial zonal wind profile of the tvpe 
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9. 


u B( 1—cos 5 y ) 


Interpreting D as the distance between the walls we 


(2.9) 


consider a profile where u=0 at both walls and has a 
value u=2B for y=D/2, i.e., in the center of the channel. 

We shall in the following use an expansion of the 
function a=a (y) in the form: 


(2.10) 


N 
aly)= > > a, sin (n\y) 
n=1 


where \=7/D. The expression (2.10) for the meridional! 
dependence of the perturbation stream function is chosen 
in such a way that y¥’ always disappears at the northern 
and southern boundaries of the region. The meridional 
velocity component r=v’ will therefore also vanish at the 
boundaries, and the boundary condition corresponding to 
rigid walls is therefore automatically satisfied. 

The expression for @ (y) is inserted in equation (2.8 
and we obtain after reduction a set of linear homogeneous 
equations which can be used to determine the possible 
The 


resulting set of equations may be written in the form 


values of ¢ and the ratios between the values of a,. 


= (u? +- (n?—4n)\?) Ba, -»+ ((u? +-n?A?) (e— B) +B Peja, 


ae 
+5 (u?+ 


(n?+4n)d?) Ba, ..=0. 2.11 


We notice in the set (2.11) that it divides naturally 
into two sets corresponding to even and odd values of » 
It should further be mentioned that we must put a finite 
value on the maximum number N which x may obtain in 
order to determine the We shall in this 
study only be concerned with small and odd values of 


n and N. 


Kliasen |7] but mostly in cases where the 8-effect, and 


values of ec. 
Larger values of N have been considered by 


naturally also the effect of divergence, was excluded by 
assumption. 

Restricting ourselves to the most simple case N=3 we 
obtain the following set of equations from (2.11) 
5 


> 


(+e Te et a B(3y?—?*) } a, 


ls 


T = {p~ + 5x" Ba 0) 


(u?— 3A?) Ba, + { (9A2+ p? 4 ye 


+B 3(9A* +-y?) -a;=0 J 





very similar, but not identical to 


also be obtained by applying the 


A set of equations 
(2.12) may naturally 
same technique as Charney [4] in deriving the stability 
criterion for a two-parameter model with horizontal as 
well as vertical shear. Instead of the series (2.10) one 


y 


obtains the equations corresponding to (2.12) by applying 
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L x10"? km. 
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Figure 1.—The imaginary part, c;, of the solution to equation 
(2.13) as a function of the zonal wavelength for a maximum 
meridional wavelength of 4000 km., B=30 m. sec.-! Curve (1): 
B=0, g@=0; (2): B=Bys, @=0; (3): B=, @=0.5X 10-8 m.-?; 
(4): B=By, @=1.0X10-" m.-?; (5): B=Bys, G=1.5X10-" m.-? 


the linearized equation (2.8) in the center of the channel 
and midway between the center and the wall. Merid- 
ional derivatives are replaced by finite differences. It is 
a matter of convenience whether one or the other pro- 
cedure is used. 

The possible values of ¢ are now determined from the 
condition that the determinant of the system (2.12) has 
to vanish in order to have non-trivial solutions. This 
condition leads to the equation: 


(A? +-y?+-q?) (9A? 4 p+ ¢q’)c? +[28(5d? +p? + ¢’) 


—1/2B (5u'+9n'+46y?d?+-5y?¢’ 
+17d*q’) je +[8?—1/28B(5y?+ 17%’) 


+1/4B?(5u'—3!+50u°A2)]=0 (2.13) 

The frequency equation (2.13) has been solved nu- 
merically for a number of cases. These cases fall into 
different series. For each value of the maximum, merid- 
ional wavelength (22/A) it was decided to vary the zonal 
wavelength (.=22/u) in each series keeping q constant. 
The effect of 8 and of divergence was tested by keeping 
8—0 and g=0 in one series of computations. The 
results of these computations are given in figures 1, 2, 
and 3. 

The maximum, meridional wavelength in the com- 
putations illustrated in figure 1 was 4,000 km. and the max- 
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Figure 2.—The imaginary part, c;, of the solution to equation 
(2.13) as a function of the zonal wavelength for a maximum 
meridional wave length of 6,000 km., B=30 m. sec! The 


curves (1), (2), (3), and (4) as in figure 1. 
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Figure 3.—The imaginary part, c;, of the solution to equation 
(2.13) as a function of the meridional wavelength for a fixed 
zonal wavelength of 5,000 km., B=30 m. sec.-' The curves (1), 
(2), and (3) as in figure 1. 


imum zonal speed in the initial profile 60 m. see.~' (B 
30 m. sec.~'). This case corresponds to a rather narrow 


intense jetstream. The abscissa in the figure is the zonal 
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wavelength in the unit of thousands of km., while the ord- 
inate is the imaginary part of the wave speed: c,;. The fig- 
ure contains five curves corresponding to the cases: (1) 
B=0, q=0; (2) B=Bys, gq=9; (3) B=Bus, ga=0.5 X10-” m.~?; 
(4) B=By, g= 1.0 10~-” m.~?; and (5) B=By, q=1.5X10-” 
i 

It is seen by comparing the two first cases that the p- 
effect has a stabilizing influence on the perturbations in the 
sense that the magnitude of ¢, is smaller for 8#0 than in 
the case 8=0 and also that the band of zonal wavelengths 
for which instability occurs is more narrow for 60. 
By comparing the first two cases with the last three it is 
furthermore seen that the effects of g, which measures 
the intensity of the divergence implied by the model, are 
to further decrease the instability of the waves in the 
model. The divergent, one-parameter model is therefore 
more stable than the non-divergent model. 

The results are verified by figure 2, which contains 
similar calculations with the maximum, meridional wave- 
length at 6,000 km. The important new information 
contained in this figure is that the waves are more stable 
when the meridional wavelength is increased. For D= 
6,000 km. it is seen that all wavelengths are stable in the 
case of B=6y and gq=1.510~-"m.~? and only the waves 
in a narrow band around 4,000 km. are unstable for 
B=, and q=1.0*10""m.~*. Figure 3 contains curves 
giving the imaginary root as a function of the meridienal 
wavelength for a fixed value of the zonal wavelength, 
L=5,000 km. 

The information contained in figures 1, 2 and 3 may 
also be expressed in terms of the time, 7’, it would take to 
double the amplitude of the unstable perturbation. This 
time can be determined from the formula 


_In2 
He; 


T (2.14) 
where c, is the positive, imaginary part of the wave speed. 
If the wavelength L=2z/y is expressed in the unit 1,000 
km., ¢; in m. sec.~', and 7 in days we obtain 
mn L 
T= 1.28 —- 


t 


(2.15) 


The curves in figures 1, 2 and 3 indicate doubling times 
which may be as short as 0.5 day, but more typically are 
of the order 1-2 days. 

The general conclusion from these calculations is there- 
fore that unstable waves are possible even when the B-effect 
and the divergence effects are included in the analysis. 
Both of these effects have a stabilizing influence on the 
waves in the sense that the waves become more stable for 
larger values of 8 and q as they do for larger values of the 
maximum, meridional wavelength. 

The general validity of the conclusions which have been 
reached in this section may be questioned because only the 
terms corresponding to n=1 and n=3 in (2.10) have been 
included in the frequency equation (2.13). If we also 
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included the next odd value, n=5, the result would be 
that (2.13) would be replaced by a cubic equation, while 
the corresponding equation would be of the fourth degree 
if also n=7 were included. A priori the possibility exists 
that several amplifying modes will be present in the solu- 
tion, and that the growth rate of these modes is different 
from the one found in this section. Such conditions are 
found in the theory of baroclinic stability (Charney [4], 
Murray [11]). In the case of barotropic stability, under 
investigation here, there exists a limited amount of infor- 
mation regarding the higher modes in the paper by Eliasen 
[7], who has found the roots in the frequency equation up 
to and including n=8. In the case of a symmetric, 
harmonic velocity profile in the zonal current he finds for 
n<8 that only one unstable mode exists, and that the 
growth rate is the same for 3 <n <6, but somewhat small- 
er for n=7 and n=8. It is therefore likely that the di- 
agrams in figures 1, 2, and 3 contain the essential features 
of the stability of the quasi-barotropic flow, although the 
problem needs further investigation. 


3. SECOND-ORDER MOMENTUM CHANGES 


It is of interest to see what we can say from the solution 
of this particular perturbation problem about the second- 
order effects of the perturbations on the mean flow. 
Such second-order effects on the mean flow will give the 
initial changes, but after a while the mean flow will have 
changed so much that the solution of the perturbation 
problem no longer applies, which means that the mean 
flow now starts to alter the structure of the perturbations 
in a way radically different from our solution. This 
particular problem will be treated by numerical calcula- 
tions in the later section. 

The second-order changes in the mean flow can be found 
from the first equation of motion in a form applicable to the 
model : 


duu , Ouv_ _—Oe 


ou 
ot TS i 


or dy of ated 


+fr. 


Defining now 


1 L 
ae i) udz 


and writing 


=() 
or 


u=U+u’, v=v’, since T= 


we obtain from (2.14) 


Ou ou’ v’ 


ot —sb 


The initial changes in the zonal wind can, therefore, 
be computed from the convergence of the momentum 
transport of the perturbations. We must next evaluate 
this quantity from the solution of the perturbation 
problem. The general solution for the perturbatio1 
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cam function can be written (see equations (2.7) and 
10)): 


,y, t)={a,"? sin Ay+a,;" sin 3dy} e#7~ 4? 


(2) 


+ {a, sin Ay+a,;™ sin 3ry}e*%-%" = (3.5) 
where ¢; and ¢, are the two solutions of the frequency 
equation (2.13) and a, a3", a, and a3” are the four 
(unknown) amplitudes of the wave components in the 
solution. The amplitudes (@) will, in general, be complex 
numbers if the wave speed ¢ is complex. Next let the 
initial perturbation stream function be given by 

Yo’ (x, y) = (A, sin Ay+A; sin 3dy)e™. (3.6) 
We have then the following four equations to determine 
the amplitudes. Two equations are determined by the 
condition that the expression (3.5) has to agree with 
(3.6) for t=0, giving: 


(4 2) 

a) . +a,$ ax As 
( ! 2 

a3 D + a! ‘= As, 


The remaining two equations are determined by one of 
the two equations (2.12), which has to be satisfied for 


(3.7) 


c=c¢, and c=. Using the first of the equations (2.12) 


we obtain 


F(e,)on +Ga"! + 
F(e2)a,° +Ga, =0 


where F(c) and G are defined by the expressions 


F(e) =(1 +X +5) c+ 5-5 B (3-%) (3.9) 


: (1 +5 *) B. 
mM 


The system (3.7), (3.8) becomes eight linear, inhomo- 
geneous equations when the amplitudes (@) and wave 
velocity (c) are written as complex numbers and the equa- 
tions are separated into real and imaginary parts. Let 
the wave speeds be written: 


and 


(3.10) 


¢=R-18, (3.11) 


¢=R-iS 
and the amplitudes: 
a, =X,4+01X, |) 
a, =Y,+1Y, 


ae (3.12) 
a;°? =Z,+ iZ, 





a;®=V,+1V;2 J 


It is obvious that the eight quantities: X), X2, Yi, Yo, 21, 
Z,, Vy, and V, can be determined from (3.7) and (3.8) 
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when F and S have been determined from the solution of 
(2.13). We may therefore consider the real and imaginary 
parts of the amplitudes of the wave components as known. 
The final expression for the perturbation stream function 
may then be written: 


VW’ (x,y, t) =e" { (X, sin Ay+-Z, sin 3dy) cos p(r— Rt) 
— (X, sin \¥y+Z, sin 3dy) sin p(x—Rt) } 
+e-S'/(Y, sin Ay+V, sin 3dy) cos u(z—Rt) 
—(¥, sin \y+ V2 sin 3Ay) sin w(x—R2) } 
(3.13) 


From the knowledge of the expression for the perturba- 
tion stream function it is now a straightforward (although 
tedious) computation to evaluate the convergence of the 
meridional transport of momentum to be used in (3.4). 
We find: 


Ou 


yi te~*S'(Y V.—Y.V)) 


=2yud?|¢ 2ySt ( X,Z2— Ask ) 


T | ViX.— V.X 14 Zi o- Z.Y ;)] -M(y) (3.14) 
where 

M(y)=cos 4y—cos 2dy. (3.15) 
It can be shown that for our case we will always have 


—_— (X,2.—X2Z, ) = ( Y.V; — Y2V;) — N* (3.16) 
and 
(3.17) 


ViX2—V2Xi+Z,Y2—Z2Y :=0. 
Using this information we may write (3.14) in the form: 
ou 


yeaa sinh (2uSt) -M(y) 


(3.18) 


or by integration 


2\7* 


“ ih N 
u(y,t)=U(y,0)—2 x [cosh (2uSt)—1]-M(y). (3.19) 


In the cases which have been treated numerically it 
has invariably shown up that N* is a positive quantity. 
This result is in agreement with remarks made by Kuo 
[8], on the basis of results obtained by Tollmien [16], and 
we shall show later that we can obtain this result (NV*>0) 
from physical considerations of the energy conversions 
between eddy kinetic and mean kinetic energy. 

Taking this result for granted at the moment we can 
then find qualitatively what will happen to the profile of 
the zona! wind. The function M(y) is plotted in figure 4 
as a function of y/D. It has a positive maximum in the 
center of the channel and negative minima on both sides. 
Since furthermore [cosh (2uSt)—1] is positive for ¢>0 
it follows that the last term on the right side in (3.19) 
has the opposite sign of M(y). We can therefore expect 
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Figure 4.—The function M(y) as defined by equation (3.15). 
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Figure 5.—The profile of the zonal wind initially and at 2, 3, and 
3.5 days computed from equation (3.19) with B=8,, g@=O0, 
zonal wavelength 5,000 km., width of the channel 3,000 km., 


B=30 m. sec.-!, A; corresponds to 5 m. sec.~!, and A;=0. 


that the initial maximum in u(y,0) will be reduced during 
the first time period after an integration has started 
while maxima will tend to appear to the north and south 
of the middle of the channel as time goes on. These 
results are illustrated in figures 5 and 6, where we, in 
figure 5, have computed the distribution of u(y,t) for 
t=0, 2, 3, and 3.5 days from equation (3.19). Figure 5 was 
computed for non-divergent flow (¢=0), B=6,, a width 
of the channel equal to 6,000 km., and a zonal wave- 
length of 5,000 km. The initial maximum wind speed 
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Figure 6.—The profile of the zonal wind initially and at 0.5, 1.0, 


and 1.5 days computed from equation (3.19) with parameters as in 
figure 5 except A; corresponding to 5 m. sec.~! 


was 60 m. sec.~! and A; was chosen corresponding to an 
initial meridional wind of 5 m. sec.~', while A,;=0. 
Figure 5 should be compared with figure 6 where the 
parameters have the same values except that A; also 
corresponded to 5 m. the meridional wind 
component. The distribution of u(y,t) is only given for 
t=0, 0.5, 1, and 1.5 day in figure 6. It is seen that the 
presence of the waves on the smaller meridional scale, 
initially, causes the splitting of the jet stream into the 
two jet maxima to take place at a faster rate. 

Figures 5 and 6 are very similar to those constructed by 
Thompson [15] from his theory of large-scale turbulence 
in barotropic flow as they, indeed, should be because the 
The changes predicted in 


see.-! in 


physical model is the same. 
the zonal profile in our computation are however, due to 
the presence of unstable barotropic waves, while Thomp- 
son’s theory is based upon the interaction between eddies 
and the mean flow expressed partly in quantities which 
can be derived from the mean flow and partly in certain 
statistical properties of the eddies. It will be noted also 
that the time scale of the changes in the mean zonal flow 
is the same in the two computations. It takes 2 to 3 days 
to develop two distinct maxima in the profile of the zonal 
wind from the single maximum present initially in the 
center of the channel. 

From our calculation we may therefore conclude that 
if we have barotropic instability in a flow characterized by 
a single maximum in the zonal wind, it is a dynamical 
consequence caused by the convergence of the meridional 
transport of momentum by the eddies that the zonal pro- 
file will change to a profile with double maximum. This 
is done by momentum being transported away from the 
center of the channel in both directions by the eddivcs, 
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which in turn means that the unstable barotropic dis- 
turbances have a tilt of the trough and ridge lines, which 
govs from southeast toward northwest south of the center 
and from southwest toward northeast north of the center 
This tilt is opposite to the one generally 
This fact points to the 


of the channel. 
observed in the atmosphere. 
importance of barolinic processes in the real atmosphere, 
but the experiences in barotropic prediction show that 
processes similar to those described in the present paper 


definitely occur. 
4. AN EXAMPLE OF BAROTROPIC DEVELOPMENT 


Figure 7 shows a section of the 500-mb. analysis from 
May 15, 1961 at 1200 cmt. The interesting feature on 
this section of the hemispheric analysis is the fairly straight 
jet stream from the eastern part of Canada over New- 
foundland toward the southeast over the Atlantic Ocean. 
Notice in particular that the isotherms are almost parallel 
to the contours over a large section of the jet stream. The 
isotachs for the map (not reproduced) show a maximum 
of about 30 m. sec.~' over Newfoundland with a drop to 
less than 10 m. sec.~' on both sides in a distance of about 
1,000 km. 

Figures 8, 9, and 10 show analyses for May 16, 17, and 
18 also for 1200 amr. Between May 15 and 16 a develop- 
ment took place along the jet stream over Newfoundland. 
The center of the circulation was located approximately 
at 45° N., 35° W. on May 16, 1200 amr. Notice again on 
figure 8 that we find no regions of strong temperature 
advection. During the following day the number of 
closed contours around the center increased, while the 
Low drifted slowly toward the southeast with a position 
at 43° N., 30° W. on May 17, 1200 cmt. The Low re- 
mained almost stationary during the next 24-hour period 
with a slight tendency toward filling. 

The forecasts made with the 500-mb. operational model 
from the initial data at May 15, 1961 at 1200 Gar are 
shown in figures 11, 12, and 13. The 24-hour forecast 
(fig. 11) shows that the barotropic model forecasted the 
development of a closed circulation. The center located 
at approximately 50° N., 35° W., is too far to the north. 
However, the barotropic development continued in the 
forecast between 24 and 48 hours and the circulation is 
centered at 45° N., 29° W. in the 48-hour forecast. The 
72-hour forecast shows a further increase in the number 
of closed contours, but the center has now moved too far 
toward the east. 

The initial jet stream has been divided into two branches 
on the observed maps and on the forecasts. Toward the 
end of the period we find a rather broad current over the 
eastern part of North America. This current divides into 
two branches, one going toward the northeast over the 
southern part of Greenland, the other bending toward the 
southeast and continuing on the southern side of the Low, 
which was developed in the Atlantic. 

The development of the flow pattern in the example 
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500-mb. analysis for May 15, 1961, 1200 amr. Solid 
lines are contours and dashed lines isotherms 


FIGURE 7. 
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Figure 8.—500-mb. analysis for May 16, 1961, 1200 Gmr. 


shows therefore a great similarity to those predicted by 
the linear theory. This type of barotropic development 
is observed frequently in barotropic forecasting and results 
However, it is not the 


usually in cold, cut-off evclones. 
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Figure 9.—500-mb. analysis for May 17, 1961, 1200 Gur. 














Figure 10.—500-mb. analysis for May 18, 1961, 1200 Gur. 


most frequent type of development in the troposphere. 
The baroclinic developments associated with well devel- 
oped divergence fields are the more common. 

It is interesting to note that Murray [11] in a recent 
investigation of baroclinic stability of stratospheric flow, 


-24-hour one-parameter forecast starting from May 
15, 1961, 1200 Gr. 


Figure 11. 


especially the breakdown of the polar night jet located 
close to the 25-mb. surface, indicates that the polar night 
jet is baroclinically and inertially stable. 
therefore, that barotropic instability of the type discussed 
here is responsible for the breakdown of the polar night 
jet. This result is contradicted by the conclusions reached 
by Boville [2] who states that “baroclinic wave develop- 
ment, fulfilling instability criteria and energy conversions 
of the Fleagle type appear to have been identified down 
The difference in the results of 
the two (independent) studies can probably be ascribed 
to the different models which were adopted, but it is 
beyond the scope of this paper to resolve the question of 
the stability of the polar night jet. 


He suggests, 


to wave number four.” 


5. ENERGY CONVERSIONS 


We shall next turn our attention to the energy conver- 
sion between mean zonal kinetic energy and eddy kinetic 
energy, the only energy conversion which is possible in the 
The prediction equa- 


model, if g=0; 1.e., no divergence. 


tion is in this case 


0¢/dt+V-Ve+6rv=0. (5.1) 

When we divide each quantity into its zonal mean 
value and deviation from the zonal mean value and intro- 
duce these expressions into (5.1) we obtain the following 
equation for the average flow: 


d¢/dt +V7-Ve"=0 


(5.2) 
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FicgurRE 12.—48-hour one-parameter forecast starting from May 


15, 1961, 1200 aur, 


This equation may also be written in the form: 


O¢/dt=d2(w'r’ )/dy? 


because 


Vv ¢ ‘ ) oy 


but due to the non-divergence of the eddy-flow; i.e., 


Ou’ . ov’ 
=) 
Or + Oy 


we obtain 
ra) Y Ou’ ) 


- _ dtu’ o (; 
rey) Oy 


, Ov’ ) 
oy (OY oy 


The last term is, however, zero, because 


do /_, ow’ do /_, ou 
ay ( ay ) ay (* 4. ” 


The mean kinetic energy per unit mass and unit area 
may in this model be written in the form: 


oe iflien ..1 77 
k==— > u dS ~oD |, (5) au 


O Js 


from which it follows that 


dK 1 (“2 n) (= 


l "D> _ or ” 
i iT: 2 ; 5.9 
dt D Jo Ov oy Ot dy D |, y dt dy, (5. ) 


Figure 13.—72-hour one-parameter forecast starting from May 15, 


1961, 1200 Gar. 


—du/dt=0 at both boundaries at all 
It is seen from (5.9) that we obtain the change in 
kinetic energy by multiplying (5.3) by ¥ and then inte- 
grating over the width of the channel. 


because 0° /OyOdt 
times. 


This procedure 
leads to 


*D = Ou’ r’ 
u 


an 


dK 1 (?- o*u'e’ : 
—v. | dy. (5.10) 


dt v oy” 


The change in time of the eddy kinetic energy may be 
found from a similar method which results in 
dK’ 1 [ox Our’ ] 5.11) 
—_ t ay. ‘ 
dt DyJo Oy sd | 


It is then natural to define the following expression as 
the energy conversion between eddy kinetic energy and 
mean kinetic energy: 


K’-K 


(5 12) 


The general expression for energy conversion in isobaric, 
non-divergent flow shows that if the zonal wind is posi- 
tively correlated with the convergence of momentum 
transport there is a conversion from eddy kinetic energy 
to the mean flow kinetic energy, while negative correla- 
tions between the two quantities give the opposite 
conversion. 


The energy conversion (5.12) will now be evaluated for 
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the amplifying barotropic wave according to the analysis 
made earlier in this paper. Restricting ourselves to this 
part of the solution we find that the convergence of 
momentum transport is given by (3.14) 


Ou’v’ 


Oy 


Saw (X os — eA er". My). (5.13) 


When we substitute this expression together with the 
expression (2.9) for 7 in (5.12) we find after evaluation of 


the integral 


K’-K} =p Be" (X,Z,—X,Z,)=—un Be*'N*. (5.14) 
The last expression gives us the sign of the quantity N * 
which was introduced in (3.16). The eddy kinetic energy 
of the unstable barotropic wave can only increase as a 
result of the conversion of energy from the kinetic energy 
of the mean flow. The energy conversion {K’. KR} is 
therefore negative, which means that N*= —(YX,Z, 
— X,Z,) must be positive. The sign of N* was of vital 
importance for the discussion of the changes in the zonal 
wind profile carried out earlier in the paper. 

It should be noted that the expression (5.14), on which 
we have based the conclusion that N* is positive, is 
evaluated without regard to the second-order changes in 
the zonal momentum. The expression (5.14), applies 
therefore only for a short interval of time initially, a fact 
which does not destroy the argument. The importance 
of the higher-order changes in the energy conversion and 
the question of when the non-linear effects become signifi- 
cant will be investigated at the end of section 6. 

It is obvious that the linearized treatment of the baro- 
tropic problem can describe the development of the pro- 
file of the zonal wind and perturbations only up to the 
time when there is considerable interaction between the 
two parts of the motion. This is seen directly in the 
expressions which we have derived for the change in the 
zonal wind (3.18) and the energy conversion (5.14) which 
will predict changes which tend toward infinity with 
increasing time. 

From energy considerations it is equally obvious that 
the source of energy for the perturbation, i.e., the zonal 
available kinetic energy, is finite, and that the growth of 
the pertubation therefore must stop after a while, and the 
The non-linear 


interactions which are necessary for such processes can 


energy conversion must change sign. 
only be incorporated by integrating a more complete form 
of the prediction equation. Lorenz [9] has recently shown a 
simplification of the non-linear vorticity which enables 
us to retain the non-linear interaction by restricting the 
fields to a few wave components. An integration of this 


tvpe will be described in the next section. 


6. EXTENDED TIME INTEGRATIONS OF BAROTROPIC 
FLOW 


The the 
equation was carried out for two reasons. It 


extended time integration of barotropic 


was de- 
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sirable to investigate the validity of the results obtained 
from the linear analysis, and the changes in the zonal 
profile, the momentum transport, and the energy conver- 
sions could be investigated as a function of time in a more 
general case. The integration can be considered as an 
extension of Lorenz's [9] original integration, which was 
made with equations which were simplified to a maximum 
extent. Guided by the results which we have obtained 
during the linear analysis we shall formulate a system 
which is general enough to allow double maxima in the 
zonal wind profile, but simple enough to give a very 
efficient computation. 

The form of the expression for the mean zonal wind 
(3.19) in the linear case suggests that an expression of the 
form: 


u(y, t)= B(t)+C(t)— Bit) cos (2Ay)—C(t) cos 4dry (6.1) 
may be used with advantage in the numerical experiment. 
The family of curves described by (6.1) for different val- 
ues of B and C can vary between strong westerlies in the 
center of the channel with easterlies near the walls and, 
on the other extreme, strong westerlies near the walls 
with easterlies in the center of the channel. The pre- 
dicted values of B and C from given initial values will 
therefore show the changes in certain types of zonal wind 
profiles predicted from the barotropic equations. The 
stream function corresponding to the mean zonal wind 
(6.1) is: 


a — eT 
v(y, t)=D( B+0C) [1 —B +5, sin (22y) +35 sin (4Ay) 


» 


(6.2) 


where D is the width of the channel and the constant of 


integration has been determined in such a way that 


y (D, t)=0. 
The solution (3.13) 
tion suggests, on the other hand, a prescribed expression 


for the perturbation stream fune- 


for this quantity of the following form: 


t) 


E(t : 
sin 


I 
‘: F, (ft 
I 


= ; Ey ; 
W'lz, y, t)= sin Ay sin kr+ P 3Ay sin kr 
F(t) 


k SAY COs ker 


Pi. : 
sin Ay cos kr+ -sin 
which allows one wave number (*) in the zonal direction 
and two wave numbers (A and 3A) in the meridional direc- 
tion. The complete stream function is then given by the 
expression 


Vy, t) +’ (4, y, t) (6.4) 


Y(r, y, t) 
Now (6.4) can be substituted in the vorticity equation 


oy 
or 


ovy _ . 
py) = J(V*y,~)—B 


in which B=df/dy is considered to be a constant, and we 
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cai by the technique developed by Lorenz [9] transform Q el R B (6.6) 
(6.5) into six forecast equations for the six amplitudes: ke x 

B.C. E,, Es, F;, and F;, which are functions of time only. 

With the notations: the six forecast equations for the amplitudes are: 


7 ak EF s— EsF 


dC 
dt — 


dE, ,T np Q-3 , 5Q+ _ 7Q- | ie 
dt | CR sG yy BP 2(Q-4 ; BF; 2(Q4 y OF aaa | 


dE; , — 3Q-1 , 160-1 vp R 
dt =A [ de Fs+3 0904 1) BF, 2(90+1) Chi— 904 if | 


dF, Trp 0-3 pr_ 50+! pp 70-1 >, _ FR pp 
dt —4| cB, xQ+1) o* 2(Q+1) BA 2(Q+ —_ i ® | 


dF; i i Bind a, Fat 
no | rO)Estx6g41 BE. J9Q41) Ei 9Q41 Es 


—2kQ| LE, F3— E3F;] 





We notice from the first two expressions in (6.7) that tion was carried out as a check on the derivation of the 
B+ C=M,=const. (6.8) System (6.7) 

ahs Let us further define the mean square of the vorticity 
This result could have been derived a priori because it py the formula , 
simply expresses the fact that total zonal momentum of =~ 
the system is conserved. 

It has been shown by Lorenz [9] that the simplified 
barotropic equations (6.7) conserve the kinetic energy and 
the square of the vorticity when these quantities are Substituting again from (6.4) in (6.11) we obtain the 
integrated over the region. It should be noted that the following expression for V: 


kinetic energy and the square of the vorticity are evaluated 


"L 
V=r5| y (Vy) *dr dy. (6.11) 


=e 1 (\?+k*)? 
as the sum of the contributions from the components ; 
which are included in the forecasts. 

The kinetic energy per unit area is in our case evaluated (6.12) 
from the expression. 


Wvder dy. 


(6.9) served quantities, A and V’, can be computed in each time 


1 1 " ” Again we can of course obtain a direct proof that V as 


*"L (°D ] 
Kk L 3 |; , Vy Vidr dy= 3 TD 


é 0/0 expressed in (6.12) is indeed conserved. The two con- 
JI1Oe/0 e 


step during a numerical integration of the system (6.7). 
We introduce the expressions (6.4), (6.3), and (6.2) in’ Jf they are conserved during the integration, or at least 
(6.9) and obtain after evaluation of the integrals show no trend for systematic increase, decrease, or violent 


+k? 
“ie 

19\2--K? tegration to be described later in this section. 
8 Kk? - (E+ Fs’). (6.10) It is one of the purposes of the numerical integration to 
investigate the energy conversion between mean kinetic 


fluctuations, we have a good check on the behavior of the 
-(E74+F/? numerical integration especially with respect to truncation 


K + +-(’ +3 C44 


errors. Such a check was introduced in the numerical in- 


A direct proof that the kinetic energy is conserved can be energy and eddy kinetic energy during an extended period 
obtained by differentiating (6.10) with respect to time of integration. We need therefore to evaluate the formula 


| substituting from the system (6.7). This computa- for this energy conversion in our simplified model. This is 
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accomplished by a substitution of (6.4) in (5.12) from 
which we obtain 


K’-K 3 (B (6.13) 


I ()- (FE, F;—E;F)). 
The energy conversion (6.13) was computed in each time 
step during the numerical integration. 

We shall next return to the fact that the total momen- 
tum of the fluid is conserved as shown in (6.8). This con- 
servative property puts a constraint on the amount of 
zonal kinetic energy which can be released to the perturba- 
tions, as shown by Platzman and Baer [13]. We may, in 
other words, consider the total zonal kinetic energy as 
composed of two parts: the available zonal kinetic energy 
and the unavailable energy. It is easy to evaluate the two 
quantities in the simple model used in our experiments. 
The zonal part of the kinetic energy is, as can be seen 
from (6.10), 

foi aie > - 
E,== (B+ ¢ +4 (BP+C"), 


» 


(6.14) 


We can combine (6.14) with (6.8) and rewrite (6.14) in 
the form: 


E=% Mi'+5 B(B—M,) (6.15) 


from which it is seen that 


(6.16) 


6E.=( BH; M,) 6B. 


i. is therefore at a minimum if B=% M, in which case we 
also have (=%M). In our case there exists therefore a 
particular profile of the zonal wind, described by (6.1), for 
which the zonal kinetic energy is at a minimum, namely 
The minimum value of the zonal 


WM, 


the one for which B=C. 
kinetic energy is evaluated from (6.14) with B=C 


and is 


(6.17) 


This is the unavailable part of the zonal kinetic energy. 
Only 3/8 of the initial momentum squared is therefore 
available for conversion to eddy kinetic energy. Sub- 
tracting (6.17) from (6.15) we can evaluate the available 
zonal kinetic energy. It becomes 


l 


EE. ag = pe a BC: ae. (B—C)’. (6.18) 
ba 2 S 


] 
> 
zonal, kinetic energy. It should further be mentioned 
that very large values of zonal, available energy are ob- 
tained in the case that !M,=0, which means that B= —C, 


; - : Dania ” 
in which case F,,¢.=5B’. We shall 


We notice here that if B=C= 5M, we get no available, 


illustrate later by 
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Figure 14.—Case 1. The zonal] velocity as a function of time and 
meridional distance over one period. Initial parameters: B,= 
30 m. sec.~!, C, 15 m. sec.-!, Eyy=25 m. sec. Ey=F;, 


F35=0, B=By;, width 3,000 km., zonal wavelength 5,000 km. 


the results of the actual calculations that large 
changes in the zonal flow take place in this case. 
Numerical integrations of the system (6.7) have been 
In the following we 
The 
integration was performed using central finite differences 
in time except initially, where one uncentered step was 
taken. The time step has been 1 hour in all the calcula- 
tions. We shall treat four different cases. The calcula- 
tions were in all cases carried out to 24 days (576 time 
steps), but since all the computations turned out to be 
periodical we shall usually illustrate only a shorter time 
The finite difference system in time was tested 


very 


carried out in a number of cases. 
shall describe some of the more interesting results. 


period. 
by computing the kinetic energy per unit area and the 
square of the vorticity per unit area. These quantities 
showed only small fluctuations during the complete inte- 
gration, indicating that the finite difference system used 
in these calculations is good enough even for extended 
integrations with the present system. A short descrip- 
tion of the results of four integrations follows. 

Case 1: As the first example we have selected a case 
where the initial flow pattern is barotropically unstable. 
The zonal wavelength was taken to be 5,000 km. and the 
width of the channel 3,000 km., corresponding to a maxi- 
mum meridional wavelength of 6,000 km. The quantities 
B and C were selected initially at 30 m. sec.~' and —15 


~' respectively, giving a zonal wind profile with a 
-1 


m. see. 
single maximum in the center of the channel of 60 m. sec 
The profile has also weak easterlies close to the walls with 
a maximum easterly wind of about 7 m. sec.~! The 
initial disturbance was defined by putting F,=25 m. sec.”'; 
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The zonal kinetic energy as a function of 
The energy 


Figure 15.-—Case 1. 
time over one period (upper curve, right scale). 
conversion between zonal and eddy kinetic energy as a function 
of time over One period (lower curve, left scale). 


while #,=F,—F;=0 initially. The disturbance defined 


by these values has a south-north direction of the trough 


and ridge lines and has, therefore, initially no transport 


of momentum. 

The changes which take place in the zonal winds 
during the forecast are illustrated in figure 14, where 
the zonal wind profile is shown as a function of time in 
The forecast turns out to be periodic with a 
We have therefore only shown 


hours. 
period close to 44 hours. 
one period. The initial maximum in the zonal wind in 
the center of the channel breaks down and is after about 
22 hours replaced by weak easterlies of 2-3 m. sec.~! 
At the same time two maxima of westerly zonal winds 
form to the north and the south of the center of the chan- 
nel with speeds of 28-29 m. sec.~' After 44 hours we are 
back to the initial zonal wind profile. 
dicts therefore an oscillation in the zonal wind profile vary- 


The forecast pre- 


ing between a single and a double maximum in the jet 
stream. 

The which this 
described in figure 15, in which the lower curve shows 


processes determine oscillation are 
the energy conversion between zonal kinetic energy and 
kinetic During the first hours of the 
forecast the energy is converted from the zonal flow to 
the eddies. 
about 8 hours, but remains positive up to about 15 hours, 
when the sign of the conversion changes for a few hours 
The upper 
curve shows the variation with respect to time of the 
kinetic energy contained in the zonal flow. 


eddy energy. 


This conversion reaches a maximum after 


and becomes zero shortly before 22 hours. 


In agree- 
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Figure 16.—Case 2. Zonal velocity as a function of time and 


meridional distance over one period. 
(fig. 14) except C, 


Parameters as in case | 


30 m. sec.~!. 


~ 


Figure 17.—Case 2. Zonal kinetic energy as a function of time 


over one period (upper curve, right scale). The energy conver- 
sion between zonal and eddy kinetic energy as a function of time 


over one period (lower curve, left scale). 


ment with the energy conversion curve we find a decrease 
in the zonal kinetic energy up to 15 hours, when it in- 
creases slightly, and then a decrease again before it finally 
starts to increase to complete the cycle after 44 hours. 
We can compute the 
kinetic energy from the formula (6.17). 
— and C,=—15 m.sec.*! we find M, 
The minimum zonal kinetic energy is therefore F, ».;. 
140.6 
equal to the actual minimum on the curve for £,, in- 


unavailable zonal 
With B,=30 


1 


amount of 


m.sec. 15 m. see. 


m.? sec.-? It is interesting to note that this is 
dicating that all the available kinetic energy in the zonal 
flow actually is being released to the perturbations. 
It also indicates that very large amounts of kinetic energy 
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may be released to the perturbations if we specify the 
initial zonal wind profile in such a way that the momentum 
In order to test this we computed a second case. 
—('’=30 m. 
The 


this 


is small. 
(ase 2: In this case we have chosen B 
1 which makes the zonal momentum .,=0. 
unavailable kinetic energy therefore in 
case, and the total amount of kinetic energy is at the 
In this admittedly extreme case we 


sec, 


vanishes 


same time larger. 
have again a maximum of 60 m. sec.~! in the center of 
the channel in the initial, zonal wind profile with easter- 


lies to the north and south of 33 m. sec.~! (fig. 16). The 


period in the forecasts turns out to be about 54 hours in 


Very violent fluctuations take place in the 
zonal winds at the different latitudes in this case. In 
the middle of the period we find easterlies in the center 
! while the westerlies in the 


this case. 


of the channel of 55 m.see. 
northern and southern portions, at the same time, have 
increased in strength to more than 30 m.sec.~! It can 
be seen on figure 16 that the forecast passes through a 
state between 13 and 14 hours in which there is no zonal 
wind, a situation which is found again between 40 and 
41 hours. 

Figure 17 (corresponding to fig. 15 in case 1) shows the 
time variation of the zonal kinetic energy and the energy 
conversion between zonal and eddy kinetic energy. 
We find as expected that the zonal kinetic energy decreases 
to zero between 13 and 14 hours and again between 40 
and 41 hours. The energy conversion is naturally zero at 
the same time after having been positive, since all the 
kinetic energy has been converted to eddy kinetic energy. 
The secondary maximum in the zonal kinetic energy in 
the middle of the period where the energy conversion 
again is zero corresponds to the maximum intensity of the 
double jet stream seen in figure 16 at 27 hours. 

Case 3: The initial zonal wind profile for this case was 
similar to case 1 in the sense that B and C initially had 
values (By=30 1. GQ=-—15 m.sec.~), 
but the dimensions of the region were increased. The 
width of the channel was taken to be 10,000 km., equal 
to the distance between pole and equator on the earth. 


the same m.sec, 


The wavelength in the zonal direction was also taken to 
be 10,000 km., making the initial flow stable according 
The initial per- 
turbation was again prescribed by setting £,=25 m.sec.~! 
initially, while all the other amplitudes were zero initially. 

Some results from this forecast are shown in figures 18 
and 19. 
of the zonal wind in the center of the channel and at 
a distance 0.2) from the wall. The zonal wind in the 
center shows a regular variation between 60 m.sec.~! and 


to the linear theory (see figs. 2 and 3). 


The first (fig. 18) shows a record of the strength 


32.5 m.see.~'! with a period of slightly more than 6 days, 
while the variation closer to the wall is between —6 
We find 


again a tendency to divide the initial single jet maximum 


m.sec.~' and +9 m.sec.~! with the same period. 


into two maxima, because the westerly wind close to the 
wall is at a maximum at the same time as the zonal wind 
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Figure 18.—Case 3. The zonal velocity in the center of the 
channel and at a distance of 0.2 D from the wall as a function of 
' C.=15 m. sec., &, 


Parameters B,=30 m. : 
10,000 


width 10,000 km., zonal wavelength 


time. see. 
Fyy= F»=0, B=By;, 


km. 


Figure 19.—Case 3. The zonal kinetic energy as a function of time 


over one period (upper curve, right scale). The energy conver- 


sion between zonal and eddy kinetic energy as a function of 


time over one period (lower curve, left scale). 


is at a minimum in the center of the channel. The energy 
conversion between zonal and eddy kinetic energy is 
plotted in figure 19 as a function of time together with the 
time variation of the zonal kinetic energy itself over one 
period (about 148 hours). The two curves have this time 
a very regular variation, where the zonal kinetic energy 
goes to a minimum equal to the value estimated from 
equation (6.17) in the middle of the period. 

Case 4: The initial specifications for this case were 
By=30 m.sec.-', Cy=—30 m.sec.~!, the width of the 
channel 10,000 km., and the zonal wavelength 10,000 km. 
We have therefore again chosen a situation where the 
total zonal momentum vanishes. The variation of the 
zonal wind in the center of the channel and at a distance of 
0.2D from the wall is shown in figure 20 as a function of 
time over a 24-day period. The zonal wind in the center 
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Figure 20.—Case 4. The zonal velocity in the center of the 


channel and at a distance of 0.2 D from the wall as a function of 


time. Parameters as in case 3 (fig. 18) except C, -30 m. see.~!. 


varies between 60 m.sec.~' and 0 m.sec.~' with a period of 
5% days, while the zonal wind at 0.2D varies between— 33 
m.sec.~' and 0 “Notice, that no zonal wind 
exists at 2% days and at intervals of 5% days thereafter. 
This picture is consistent with the curves in figure 21 where 
we notice that the kinetic energy of the zonal flow goes to 
zero in the middle of the period. 

The characteristic period of 5—6 days which has been 
found in these integrations agrees with the typical time 
scale found by Thompson [15] in his theory of long- 
period variations in barotropic flow. A further agreement 
is found between the results obtained here and 
reported by Charney [4] from a single extended integration 
of the barotropic vorticity equation and by Baer [1]. 

Some support is also found in analysis of atmospheric 
data for periodicities of this order of magnitude. Refer- 
ence is made to the observational studies by Mintz and 
Kao [10] who found a period of 3—5 days in the convergence 
of the meridional transport of momentum at certain 
latitudes during January 1949. The study by Duggan [6] 
of the meridional convergence of momentum also shows 


m.sec, 


those 


periodicities of this order of magnitude. 

The results of the non-linear integrations can be used 
to obtain information about how soon the higher-order 
effects become important. One may, for example, obtain 
this information from curves giving the energy conversion 
from zonal kinetic energy to eddy kinetic energy as a func- 
tion of time in the non-linear computation and estimated 
from the linear theory including terms of different orders. 

Substituting the expression (3.19) for u (y, t) and (3.18) 
for —d (u’ v’)/dy in the formula (5.12) for the energy con- 
version we find first that 


= "D 
K.K’ } =4p2N* sinh (2uSt)- LJ Ty, MUydy | 


$A? : *D - 
AS - (cosh (2uSt)—1) sinh (2uSt)- E | M*y)dy |: 


(6.19) 


The first term on the right hand side of (6.19) gives the 
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The zonal kinetic energy as a function of 
The energy 


Figure 21.—Case 4. 
time over one period (upper curve, right scale). 
conversion between zonal and eddy kinetic energy as a function 
of time over one period (lower curve, left scale). 


contribution which depends upon the initial profile u 
(y, 0), while the second term is independent of the initial 
profile. Since the meridional distribution of the conver- 
gence of momentum transport, expressed by the function 
My), is a second-order effect, we may say that the first 
term in {K - K’} is of second order, while the last term 
is a fourth-order effect. 

When the expressions (2.9) for u (y, 0) and (3.15) for 
M(y) are introduced in (6.19) we obtain after evaluation 
of the integrals 

K- K’} =2u?N*B sinh (2uSt) 

Sut N* 

~ § 


sinh (2uSt) - (cosh (2QuSt)—1). (6.20) 


The following case was selected for comparison, The 
zonal wavelength was chosen to be 5,000 km., the maxi- 
mum meridional wavelength to be 6,000 km., B= 16> 
10-"m.~'sec.~!, B=30 m. sec.~!, and g=0. The curve 
without marks in figure 22 gives the energy conversion in 
the non-linear computation as a function of time. The 
to a maximum, which is 
The curve in figure 


energy conversion increases 
reached after approximately 44 hours. 
22 marked with circles is the energy conversion computed 
from (6.20) including, however, only the first term. By a 
comparison of the two curves we find that the curve 
including the second-order effects gives a good estimate of 
the energy conversion up to approximately 24 hours, 
although the estimate is slightly too high. When the 


fourth-order effects are included in the evaluation of (6.20), 
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Figure 22.—The curve without marks gives the energy conversion 
as a function of time in a non-linear calculation with the following 
maximum meridional 


1 B=30 m. see 


parameters: zonal wavelength 5,000 km., 
wavelength 6,000 km., B= 16>» 
The 


conversion computed from the linear solution using the same 


10-” m.~! see 


marked with circles gives the energy 


and ¢@ 0 curve 


parameters and including second-order effects, while the curve 
marked with squares also includes fourth-order effects 


we obtain the curve in figure 22 marked with squares. 
This curve gives a very good estimate for a little longer 
period than the second-order curve, but then rapidly 
shows an underestimate. The main result of the inclu- 
sion of the fourth-order term in (6.20) is, however, that 
this term actually predicts that the energy conversion 
will reach a maximum. The main conclusion from figure 
22 is that the higher-order effects become important after 
about 1 day. 


7. SUMMARY 


The first sections of this paper contain a stability 
investigation of the divergent, one-parameter model. 
The main conclusion is that the introduction of a diver- 
gence term into the model tends to reduce the instabilities 
present in the non-divergent model. The second-order 
effects of the disturbances on the profile of the zonal 
wind are investigated. It is found that the unstable 
barotropic disturbances tend to cause a decrease of the 
zonal winds in the center of the channel and an increase 
of the mean zonal winds to the north and south. 
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A simple non-linear barotropic model containing only 
a few wave components is investigated in the next sections 
through integrations over extended time periods. The 
changes in the profile of the zonal wind, the energy con- 
version between kinetic energy of the mean flow and the 
eddies, and the time periodicities of the flow are investi- 
gated. Fluctuations of the order of 2 days are found in 
initially unstable flow, while the large-scale stable flow 
pattern shows periods of the order of 6 days. 
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ABSTRACT 


During the period of the Nashville Community Air Pollution Study, the Weather Bureau at Nashville, Tenn., 


forecasted air pollution to be one of three levels for the following day. 
Discrimination of air pollution levels can undoubtedly be im- 


the winter season when pollution ranges were high. 


The forecast method was most reliable during 


proved through experience and further studies, and the method is general enough to be used at other locales. 


1. INTRODUCTION 


An air pollution forecast was provided daily for the 
City of Nashville, Tenn., from January 23, 1959 to Au- 
gust 1, 1959. The forecasts, which were issued by the 
U.S. Weather Bureau at Nashville, were for high, low, and 
little or practically no pollution for the following day. 
The forecast method was intentionally made simple so 
that specialized experience would not be required by the 
forecasters and also took into account the type of routine 
which The criteria 
were to some extent based on past studies and experience 
of meteorologists at the Weather Bureau Research Sta- 
tion, Cincinnati, Ohio. This was probably the first time 
that a Weather Bureau Office made daily quantitative 
forecasts of air pollution levels. 


forecasts could be made. selected 


2. METHOD 


The forecast pollution level (i.e., high, low, little or none) 
was objectively determined from forecast values of a sta- 
bility index and surface wind speed. Meteorologists at 
Berry Field, approximately 8 miles southeast of Nashville 
proper, could make these forecasts without prior knowl- 
edge or experience of air pollution meteorology. The 
low-level stability index was obtained by algebraically 
subtracting the forecast surface temperature for 0600 cst 
tomorrow from the forecast 900-mb. temperature for the 
This index value, together with the forecast 
of weather the 0600 cst surface wind speed would be 
greater or less than 7 kt., was used to forecast the inten- 
sity of the air pollution as follows: 


same time. 


Forecast air pollution 
level 


Forecast stabilily index (SI) and 
surface wind spced (WS) for 
0600 cst 


‘0 and WS<7 k 
<0 and WS<7 kt. 
or : Low pollution 
‘0 and WS>7 kt. 
0 and WS>7 kt Little or no pollution 


High pollution 


The forecast was made at approximately 1500 cst and 
was valid for the next morning. Wind speeds were fore- 
cast using prognostic positions of pressure systems and 
The 900-mb. temperature 
forecast was based on 850-mb. up-wind radiosonde obser- 


prognostic pressure gradient. 


vation information, and wind patterns and temperature 
advection as indicated on the 850-mb. chart. The latter 
was adjusted to the 900-mb. level, keeping in mind lifting 
or subsidence, as indicated. Surface temperatures were 
forecast using the 850-mb. temperature forecast and ad- 
justing it to the surface with due allowance for forecast 
cloud cover or the lack of it, as an indication of radiational 
cooling and with allowance for wind speed as an indication 
of vertical mixing at low levels. 


3. RESULTS 


For each day, the air pollution at Nashville is described 
by each of the following: 

A. Forecast Pollution Level, FPL:—The level of pollu- 
tion forecast by the Weather Bureau at Nashville using 
the previously described method. 

B. Pollution Level From the Observed 
Criteria, PL(OC):—The stability index and wind speed 
which actually occurred determined this level in the same 
way that the the FPL. 
In other words, would be 


Mete orolog ical 


forecast elements determined 
the PL(OC) and the FPL 
identical if the meteorological elements were correctly 
forecast. 

(. Measured Pollution Level, MPL:—In conjunction 
with the Nashville Community Air Pollution Study, there 
were 32 sampling stations uniformly spaced over the 


city from which continuous measurements of sulfur 
dioxide (pphm) and soiling index (Coh/1000 linear ft.) 


were available. Coh values determine by filter paper 
method, are an approximation of the soiling properties of 
the atmosphere. At seven of these stations, which were 
also uniformly spaced, more extensive sampling was done, 
and continuous 24-hour total particulate concentrations 
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OR 
NONE NONE NONE 
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Figure 1.—Total particulate compared to pollution level indicated 
by observed meteorological criteria, PL(OC), separated by seasons 
and air quality group. Horizontal line is geometric mean, ver- 

tical line gives limits of one standard deviation. 


(ug./m.*) were measured using high-volume air samplers. 
Daily means for each of the three air quality measure- 
ments were computed and used as a basis for determining 
the MPL’s. 

Comparisons of the three descriptions of “pollution 
level’ were thus possible. A measure of ability to forecast 
the meteorological criteria was obtained by a comparison 
of the FPL and PL(OC) for each day a forecast was made. 
The effectiveness of the meteorological criteria for indi- 
cating pollution level was evaluated by a comparison of 
the PL(OC) and the MPL for each day. Finally, the 
FPL and the MPL were compared to determine the merit 
of the program. The results of the three comparisons 
are discussed in the following subsections. 


EVALUATION OF FORECAST ABILITY: FPL VS. PL(OC) 


Table 1, a contingency table of FPL vs. PL(OC), 
reveals that the forecasts were right more often than 
wrong with the exception of the “little or none’ forecasts. 

A statistical test (x?) indicates that such a distribution 
has less than 0.1 percent probability of occurring by 


chance. For the entire period 45 percent of the forecasts 


TABLE 1.—Frequency of forecast air quality level (FPL) classified 
according to the air quality level from the observed meteorological 
variables, PL(OC). 


High Little or Total 
none 


High 
Low 
Little or none 


Total 


Calculated x?= 23.3 
From tables x? (.999 
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OR 
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SO, measurement compared to PL(OC), separated into 
seasons and air quality groups. 


Figure 2. 


were meteorologically correct. If the period is divided 
into seasons (table 2), then 53 percent of the forecasts 
were correct in winter, 27 percent in the spring, and 53 
percent in the summer. The low springtime figure 
probably reflects the erratic weather conditions during 
this season. For the purpose of the study, winter was 
presumed to be from January 23 to March 15, spring from 
March 16 to May 15, and summer from May 16 to 
August 1. These dates were primarily determined from 
heating degree days although an attempt was made to 
keep the number of days in each season somewhat equal. 


EVALUATION OF EFFECTIVENESS OF METEOROLOGICAL CRITERIA: 
PL(OC) vs. MPL 


Because of differences in range and distribution of air 
quality data from one season to the next, comparisons 


were made season by season. Each group of air quality 


TABLE 2. 


FPL vs. PL(OC) by seasons 


PL(OC) 
Percent 
correct 


High Little or 
none 


Winter 


High_. 

Low - 

Little or none 
Total 


Spring 


High... 

Low 

Little or none 
Total 


Summer 


High. 

Low... 

Little or none 
Total 
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Ficure 3.—Soiling index values compared to PL(OC), separated 


into seasons and air quality groups. 


data was tested for normality of distribution by season, 
and it was found that all air quality data for all seasons 
were distributed in the logarithmic normal. For each 
season the three types of data, i.e., soiling index, sulfur 
dioxide, and total particulate, were separated according 
to the PL(OC), giving three different groups of data for 
each air quality parameter. The geometric means and 
standard geometric deviations from the groups thus 
obtained were plotted on linear graph paper. In figures 
1, 2, and 3 the horizontal line is the geometric mean and 
the vertical line represents one standard geometric devia- 
tion on either side of the geometric mean (i.e., approxi- 
mately 68 percent of the data are between these limits) 
Each parameter and season should be considered sep- 
arately, but the one thing that all have in common is 
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that the means are in the proper order. That is, the 
data selected on the days which were classified as high 
pollution according to the observed criteria had higher 
means than those of the low pollution days, and so on. 
The figures also indicate very definite seasonal differences 
which vary somewhat according to the type of measure- 
ment. Coh and SO, values are high with a wide range 
during the winter months, whereas the summer values 
are too low to be considered. 

A statistical test, the analysis of variance method, was 
emploved to determine if the means differed significantly. 
The test indicates very significant differences among the 
means during the winter season for all three air quality 
varameters. Apparently the method effectively divides 
the data into three levels during this season. The SO, 
values during the spring and summer seasons are generally 
too low to be considered. Total particulate as well as 
soiling index means during the spring season differ sig- 
nificantly, but an inspection of the figures reveals that 
good discrimination exists between only two levels. For 
the summer season there are significant differences among 
the total particulate means, but not as marked as during 
the winter season. Significant differences among the 
soiling index means are indicated, but the range of the 
data is low. To summarize the three figures, it appears 
that the criteria are valid indicators of three levels of 
total particulate, soiling index, and sulfur dioxide during 
the winter season. During the spring season the criteria 
signiiicantly indicate two levels of total particulate and 
soiling index, while during the summer season they sig- 
nificantly indicate three levels of total particulate. 

A different approach was also used to test the effec- 
tiveness of the meteorological criteria. The frequency 
distributions were divided into thirds, and each day’s 
data were designated as being high, low, or little or none, 





TaBLeE 3.—PL(OC) vs. MPL by seasons 
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é | ages of agreement were computed. As expected, these 
a adien . . . . 
- J figures are highest during the winter season when weather |, 
w 30r— — ad ~—e = . 
2 conditions are sharply defined and pollution ranges are — { 
$ high. With an overall winter average of 70 percent, the i 
20r— on \9 ° ° . ° 
Ss PL(OC)’s agree with the air quality level, as determined | ; 
a ° ° ° ° a 
3 from the frequency distribution, approximately 7 out of | | 
” . rm ° 
a | 4 10 times. The averages during the other two seasons are 
T significantly lower with the highest percentages of agree- 
a —_ | ment for the high pollution and little or no pollution ‘ 
HIGH Low LITTLE =— WIG Low LITTLE a , a _* . 
x3 ss. levels. The criteria effectively indicate three levels of 
FPL pollution during the winter season and only two levels ' 
Figure 5.—SO, values compared to FPL. during the spring and summer seasons. ( 
- , | 
TaBLeE 4.—F PL vs. MPL by seasons | 
MPL | 
FPL weg./m.? SO; Coh/1000 linear ft. 
High Low Little Percent High Low Little Percent High Low Little Percent 
or none correct or none correct or none correct 
WINTER 
High . 3 1 66.6 10 1 1 83.3 s 3 1 66. 6 
Low... 3 6 3 50.0 3 5 4 41.7 3 3 6 0 
Little or none 1 1 4 66. 6 0 2 4 66. 6 0 3 3 0.0 
ye ag 60.0 . 63.3 46.7 
SPRING 
High... : 6 1 2 66. 6 4 + 1 44.4 3 5 1 33.3 
eae ; 6 y 0 60.0 6 5 4 33.3 6 4 5 26.7 
Little or none. 3 3 7 53. 8 2 5 6 46. 2 3 3 7 3.8 
a oii 59. 5 _ 40.6 7.8 
SUMMER 
High... . piacoaus 9 6 1 £6. 2 
Low # a 8 8 6 36.3 Data range too low Data range too low 
Little or none... ° , » 0 4 7 63.6 
Total. ... 49.0 
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EVALUATION OF PROGRAM: FPL vs. MPL 


similar methods were used for the evaluation of the 
forecasts in terms of the measured pollution level. Here 
again, the data from the air quality measurements were 
separated first by season and then according to FPL. 
The geometric means and standard geometric deviations 
from the distributions thus obtained were plotted (figs. 
4,5, and 6) and compared. Again it is noticed that the 
means for all seasons are in the proper order. ‘The 
forecasts divide the total particulate and sulfur dioxide 
data into three significantly different levels during the 
winter season, while the soiling index data are divided 
into only two distinguishable levels. During the spring 
season the SO, values are again too low to show discrimi- 
nation and there is only recognizable discrimination 
between two levels of soiling index and total particulate. 
Soiling index and SO, data are too low to be considered 
during the summer season. Three significant levels of 
total particulate are, however, indicated, although the 
associated large standard deviations indicate that the 
discrimination is poor. 

Table 4 is similar to table 3 except that the FPL is used 
instead of the PL(OC). The period of highest forecast 
verification occurred during the winter season with ap- 
proximately 6 out of 10 forecast pollution levels in agree- 
ment with the air quality levels as determined from the 
frequency distributions. The percentages, however, do 
indicate that the winter forecasts are more accurate for 
the highest and lowest levels as compared to the middle 
level. During the spring and summer seasons approxi- 
mately half of the forecast levels agree with the associated 
total particulate levels, and the percentages indicate a 
greater agreement for the highest and lowest levels during 
the summer season. SO, and soiling index percentages 
are somewhat lower during the warmer seasons because 
of the low range of the data. This is a consequence of the 
increased relative importance of ordinarily random varia- 
tions of source strength, wind direction, local currents, and 
other meteorological variables. 
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4. CONCLUSIONS AND RECOMMENDATIONS 


The tested meteorological criteria are shown to be good 
indicators of three air pollution levels during the winter 
season, but can discriminate between only two levels dur- 
ing the spring and summer. Of the three methods of 
measuring air quality, only total particulate maintains 
high enough values to be considered as a year-round air 
pollution indicator. Sulfur dioxide and soiling index 
measurements are very low during the summer season and 
considered alone imply the absence of community air pol- 
lution problems during that season. 

While the results are not outstanding, it is believed that 
forecasting skill would improve with further experience. 
More elaborate meteorological criteria could have been 
used which might have given better discrimination among 
the levels, but such a method would lack the desired forth- 
right manner of the one used. Other straightforward 
criteria covld be devised and tested using the same air 
quality data and associated meteorological conditions. 
Among the various weather parameters which could be 
investigated are average wind speed, wind direction, and 
the presence or absence of an upper inversion. An addi- 
tional influential parameter could be the duration of low 
wind speeds and/or interruption of these low speeds by 
periods of significant wind speed. 

The primary value of this program has been to demon- 
strate that objective forecast methods can be used by 
local Weather Bureau offices to forecast daily air pollution 
levels during those seasons of the year when air pollution 
concentrations are high and troublesome. The forecast 
method is general enough that similar methods could be 
applied to other cities after due consideration of climatol- 
ogy, topography, industrial locations, and the individual 
needs of the organizations which are using the forecasts. 
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THE OCCURRENCE OF ANOMALOUS WINDS AND THEIR SIGNIFICANCE 


M. A. ALAKA 


ABSTRACT 


Observational evidence is provided for the occurrence of anomalous winds which represent an anticyclonic rotation 


in space, and a mechanism for their development is suggested. 


The unstable nature of these winds and the role they 


play in the development of certain types of atmospheric disturbances is then discussed, and it is suggested that anoma- 


lous winds provide the dynamic mechanism for triggering hurricane formation and for the observed deepening of 


troughs downstream from intense pressure ridges. 


Finally it is noted that although the observational evidence pre- 


sented is for the occurrence of anomalous winds over small regions of the atmosphere, their development is dependent 


on large-scale processes and their effect extends beyond the area where they occur 


1. INTRODUCTION 


The gradient wind equation is a quadratic and thus has 
two solutions. One of these solutions, appropriate to anti- 
cyclonic flow, represents a clockwise rotation in space and 
is, therefore, in the opposite sense to the earth’s rotation. 
Meteorologists have traditionally given little attention to 
this solution and some even consider it as an algebraic 
Although 


there are cogent reasons to believe that winds correspond- 


accident with little or no physical significance. 


ing to this solution do not occur on a large scale, the con- 
tention, ipso facto, that this solution is of no importance is 
not justified and has probably had far-reaching effects in 
eliminating from consideration some promising avenues of 
research in connection with the development of atmos- 
pheric disturbances. 

In recent vears there has been a slow trickle of evidence, 
both theoretical and observational, that these so called 
“anomalous winds” occur more frequently than had been 
suspected. But, by and large, meteorologists are still not 
sufficiently aware of these winds nor of their importance. 
The purpose of the present study is not only to provide 
further evidence of the occurrence of anomalous winds, 
but also to demonstrate that these winds are germane to 
the development of dynamic instability in curved airflow 
and, as such, are relevant to the development of certain 
tvpes of atmospheric circulations. 


2. ELEMENTARY DYNAMICS OF ANTICYCLONIC 
MOTION 


To introduce the problem, we shall review some elemen- 
tary aspects of the dynamics of anticyclonic motion. 
The gradient wind equation may be written: 


(1) 


f=Coriolis parameter; and 6, 


were K,=1/R, 
for cyclonic motion; V 


trajectory curvature considered positive 
magnitude of the horizontal wind; 
— adp On 


pressure gradi- 


ent force. From equation (1) 


3, (lt 2) 


> 


F 


\ 4K,b, ) 


Since V must always be positive, the solution with the 
plus sign has a physical meaning only when A, is negative; 
i.e., When the flow is anticyclonic. 

We shall confine our discussion to anticyclonic flow and 
put K;=—K,. The two solutions of equation (2) reduce 
to: 

V ft 1 1 4K'b, ) 
An ar 
and 


. 4K‘), 
Ving (1ty 1p) 


From the above equations it can be seen that 


We differentiate equation (1) with respect to V 


db ort a 
"—IVK,- () 

d\ ts 

and note that db,/dV can be zero; i.e., 6, has an extreme 


value if A, is negative. This value is reached when 
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Figure 1.—Variation of gradient anticyclonic wind speed with 


pressure gradient. 


a ar 
2K, 2K; 


and represents the maximum value 


eee 
4k, 4K; 


which the pressure gradient can attain in anticyclonic 
motion, 

Figure 1 shows the variation of V with 6, in accordance 
with equation (1). It is seen that the two solutions meet 
at \'=//2K,’ when the pressure gradient is at its maximum 
value P/4Ky’. 
to twice the geostrophic wind V,, the two solutions V; and 


From this point, where the wind is equal 


\V,, corresponding respectively to equations (3) and (4) 
V7, decreases with decreasing pressure gradi- 
This 
corresponds to the conditions normally observed in the 
We shall therefore term V, the normal solu- 
V2, on 


branch out. 
ent and vanishes when the latter becomes zero. 
atmosphere. 
tion of the anticyclonic gradient wind equation. 
the other hand, increases with decreasing pressure gradient 
and at 6,=0 reaches the value V,=//K,’ and the flow 
becomes inertial. Following Gustafson [7], we shall term 
V, the anomalous solution. 

Figure 2 shows the variation of the gradient anticy- 
clonic wind with trajectory curvature, assuming a con- 


Stunt pressure gradient. It is seen that, in contrast with 
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Variation of gradient anticyclonic wind speed with 
V, denotes the geostrophie wind. 


FIGURE 2. 
trajectory curvature. 


normal anticyclonic winds which increase with increasing 
curvature, anomalous winds show a sharp decrease with 
increasing curvature. 

As mentioned in the introduction, meteorologists in 
general have attributed little importance to the anomalous 
solution on the basis that it is seldom, if ever, realized in 
the atmosphere [5] or that its occurrence is limited to small- 
scale mechanically produced vortices, or to atmospheric 
eddies produced by friction [11]. The standard arguments 
in support of this view are usually some variation of the 
following: 

a. The anomalous solution requires a clockwise rotation 
in space and, therefore, in the opposite sense to that of the 
earth. There is no known mechanism capable of pro- 
ducing such a motion on a large scale. 

b. By expanding the quantity under the radical sign 
in equations (3) and (4), we obtain 
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, f 2K, 2K,,', 3. EVIDENCE OF THE OCCURRENCE OF ANOMALOUS 
Visser | 1-(1 et ---) WINDS 
2K, : he 
b,  Kib,? In contrast with the above arguments, there are others 
fr p re (19) which indicate that both solutions for gradient anticy- 
clonic flow are satisfied. Godson [6], for instance, con- 
and siders a moving anticyclonic streamline system and ex- 
__ f [ 1+( 12K be_2Ki%.' np ) presses the gradient wind relation as follows: 
2K; i I P 4 C cos RO ‘ - 
; Pa V J Pr tees¥ sf (/ at yt) - IRV, | (12) 
, } 2. 2 2 
ann (11) 


aes FF 


V, is thus continuous when the curvature decreases in- 
definitely and the isobars are nearly straight. On the 
other hand, for indefinitely small curvature, V, becomes 
infinitely large and would require an infinite supply of 
energy, Which is not available in the atmosphere [5]. 


contour curvature, C=speed_ of 
angle from the direction of motion 


Godson 


where R, of 
streamline system, y 
of the streamline system to the wind direction. 
argues that if C cos Y<0, it is possible to obtain two values 
of V both of which are greater than 2V, and thus 
On the other 


radius 


correspond to the classically rejected root. 
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FIGURE 4. 


hand if C cos y>0, it is pe ssible to obtain values of V both 
of which are less than 2V, and thus correspend to the 
classically accepted root. He concludes that both roots of 
the gradient wind equation for anticyclonic flow occur in 
the atmosphere. 

On the observational side, indirect evidence that anom- 
alous winds occur was provided by Gustafson [7] who 
compared observed winds with geostrophic winds in a 
moving pressure ridge at 700 mb. and applied the following 
criteria which follow directly from equations (5) and (6) 


abov PS 


ani (1:33) 





(NAUTICAL MILES) 


Wind field at 35,000 ft. around hurricane Daisy on August 26, 1958 


where 8 is the deviation in direction between the observed 
and geostrophic winds, measured clockwise from due east. 
In a 3-day sequence Gustafson found consistent positive 
deviation from geostrophic directions to the west of the 
moving ridge and similarly consistent negative deviations 
to the east of it. 

More recently, Angell [3] has furnished more direct 
evidence of the occurrence of anomalous flow from an 
evaluation of transosonde flights at 300 mb. In 
instance the transosonde performed a full clockwise loop 
in 36 hours with an estimated angular velocity of —0.48 
10-* sec.-! This represents anomalous flow since the 
angular velocity of the earth at that locality (30°N.) is 
0.36 10-* sec.-' More significantly, perhaps, the 
transosonde trajectories suggested the occasional occur- 


one 
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rence of anomalous flow on the ridges of long waves in the 
westerlies. The significance of this fact will be discussed 
in section 7, below. 
the 


Research Project 


smce 


the 


U.S. 


summer of 1956, National Hurricane 
(NHRP), Weather Bureau, has 
been operating three specially instrumented airplanes to 
make detailed observations in and near hurricane cores. 
A discussion of the characteristics and properties of the 
instrumentation was given by Hilleary and Christensen 
{10} and will not be repeated here. Among the most 
successful missions flown were those in connection with 
hurricane Daisy which developed near the Bahamas on 
August 24, 1958. This hurricane had a well-defined and 
concentrated wind circulation and presented a clear-cut 
radar configuration which greatly facilitated the location 
of the storm 


core. Flight missions were made at dif- 


ferent levels on four days from August 25 to August 28 
1958. Among the elements measured were the wind, the 
Quasi- 


instantaneous values of these parameters were punched on 


’ 


temperature, and the radio and pressure altitudes. 


40 2 


Wind field at 34,200 ft. around hurricane Daisy on August 27, 


060 85 6 


(NAUTICAL MILES) 


o 2 


1958. 


cards at specified intervals ranging from 10 seconds away 
from the core, to 2 seconds in the core. The punched 
cards were then evaluated by machine processing and the 
various parameters were plotted on a coordinate system 
fixed with respect to the storm center. 

Figures 3, 4, and 5 show the wind field obtained by 
analyzing the observations made in the upper troposphere 
on August 25, 26, and 27. In making these analyses, it 
Was necessary to regard as quasi-simultaneous, observa- 
tions which, in fact, were made over a period of several 
hours. It is, however, believed that this shortcoming 
does not invalidate the results obtained. 

From the wind fields thus obtained, values of the speed 
and direction were plotted on a rectangular grid of points 
20 n.mi. apart. From these, computations of the quan- 
tity 2VK, were made on the IBM 650 computer. The 
curvature of the trajectory K, was computed from the 
following relations: 

If uw and v are the westerly and southerly components 
of the wind, and if we define 
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Figure 6.—Field of the quantity 2V/R, at 35,000 feet (pressure 
altitude, U.S. Standard), around hurricane Daisy on August 25, 


1958. Areas with anomalous winds are shaded. 


then 


dy ( udv_idu ) 


K, = a > =cosy 


1 dv du 
=("Z-") 


l dv d 
(“aa ) 


ds r ds 


(15) 


If the motion of the storm is represented by the vector ¢ 
with components ¢, and ¢, in the east and north directions, 
respectively, and if we assume steady state conditions, 


a ' 
K- V [u(W—e)-Ve—v(V—e)-Vu] 


F { ul(u—e) a+ (v—e, 
+ 
—?P | (u —Cz) ou + (rey) ot 


Figures 6, 7, and 8 represent fields of the quantity 2K, 
V=2 V/R, corresponding to the wind fields of figures 3, 


(16) 


4, and 5, respectively. Regions where this quantity is 
negative and numerically greater than the Coriolis pa- 
rameter, 


i.e., areas where the winds are anomalous, are 


shaded. The consistency of the patterns of these areas 
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Field of the quantity 2V/R, at 35,000 ft. for August 26, 
1958. 


FiGuReE 7.- 


from day to day attests to their authenticity and demon- 
strates beyond any reasonable doubt that anomalous 
will be noted, however, that the 
areas covered by these winds represent narrow strips of 


winds do occur. It 


the order of 1 to 2 degrees of latitude in width, and are, 
therefore, likely to escape detection by ordinary synoptic 
analysis. In addition, the vertical extent of these winds 
is also limited, as can be noted from the fact that they 
are not found in figure 9 which represents the field of 
21/F, a few thousand feet below that of figure 6. This 
is perhaps the reason why meteorologists in general have 
failed to realize both the reality and importance of these 
The latter stems from the fact that they repre- 
sent an important mechanism for the development of 


winds. 


atmospheric instability as we shall now proceed to show. 


4. THE CONCEPT OF INSTABILITY OF 
ATMOSPHERIC MOTION 


The concept of instability in the atmosphere appeared 
with the initial developments in the theory of atmos- 
As early as 1878, Rayleigh [13} 
investigated conditions under which small displacements 


pheric disturbances. 


at the boundary between two air streams grow into larger 
At about the same time, Helmholtz [S] 
found that, under certain conditions, unstable gravity 


disturbances. 


waves developed at the interface between two liquids of 
different density and velocity. 
nated in the polar front theory of cyclone formation. 


The above studies culmi- 





MONTHLY WEATHER REVIEW 





; v mi T 


HURRICANE DAISY 34,200 FT. 2yv 
AUGUST 27, 1958 Ry 


M2 Eee T T 7 7 T T 


(10* sec’) 


o 8 § $ 8 


DISTANCE (NAUTICAL MILES) 
ny 
° 


2. 2 2. 











Ry Se Tee 
40 6120 00 86 © 40 0 0 BD 4 6 8 100 120 


DISTANCE (NAUTICAL MILES) 


Field of the quantity 2V/R, for August 27, 1958, at 


34,200 ft. 


Ficure 8. 


Paralleling the above studies, another series of investi- 
gations, beginning with those of Helmholtz [9] and Ray- 
leigh [14] attempted to link instability in air currents with 
the quasi-horizontal distribution of the kinematic proper- 
ties of the currents. This type of instability is usually 
known as dynamic instability. 

In recent years, attention has been redirected to still 
another type of instability known as baroclinic instability. 
This type of instability, the importance of which was first 
pointed out by Margules [12], occurs by virtue of the 
vertical wind shear. 

There is little doubt that the above three types of 
instability all play an important role in the atmosphere 
and there is an obvious need for a unified theory which 
combines them and brings out their relative importance 
and their interrelations. In the absence of such a theory 
we must consider the particular type of instability which 
appears to have the most direct bearing on the problem 
at hand. 
lous winds is best demonstrated by studying the cireum- 
stances attending the development of the so-called 
dynamic instability, with especial reference to atmospheric 


In the present case, the importance of anoma- 


currents in curved motion. 
5. DYNAMIC INSTABILITY IN CURVED AIRFLOW 
In 1936, Solberg [16] investigated the conditions under 


which dynamic instability occurs in a steady symmetrical 
polar vortex which is initially in gradient equilibrium. 
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Ficure 9.—Field of the quantity 2V/R, at 15,600 ft. for August 25, 
1958. 


Solberg’s approach may be extended to any steady vortex 
in which the tangential variations of the wind are negli- 
gible in comparison with those in the radial direction. 

Let us then consider such a vortex in which the angular 
speed w is a function of the distance ? from the axis of 
rotation. We shall proceed by the standard method and 
impose a small perturbation on this steady vortex and find 
out the circumstances under which the frequency of the 
total perturbed motion becomes imaginary. 

If we disregard friction, we can express the equilibrium 
of forces acting on the steady vortex by the following 


equation, 


where V is the magnitude of the speed, 2 the earth’s rota- 
tion, @the potential temperature, g the gravity vector, and 
(p/1,000) “/* 

Let us now apply a uniform radial impluse to all the 
particles at distance 7’ from the axis so that at the time 
t=t the particles form another circle with radius R. We 
stipulate that at any instant, the perturbed particies 
immediately acquire a pressure equal to that of the points 
which they are occupying at that instant. We further 
assume that the perturbed particles, in their displacement 
over a small distance S from their equilibrium position 
conserve both their potential temperature 6’ and _ their 


T=C, 
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absolute angular momentum M/{=V’R+QR” so that 


6—0’=$ - Ve 
(18) 


ana 


M.—M‘=S-VM, 


In addition to the forces indicated in equation (17), the 
perturbed particles sustain an acceleration d’$/dt? so that 
their motion is determined by the equation 


Ps VY” 


de? R ] (19) 


From equations (17) and (19), we obtain by subtraction 


’gS V’—V"” 22 
dt? R? ai R 


Now 


(V—V’) R+(—0’)Vr=0 (20) 


V?—V"’=(V—V’) -2V—(V)’ (21) 
where 6 V=V—V’ is a perturbation quantity, the square 
of which may be neglected. Equation (20) thus becomes 


aS 2 wo 2)(V—V’) |R+(9—0’)Vr=0 (22) 
dt? R , 


We now express the above equation in terms of the con- 
servative quantity 17, and make use of the relations in 


(18): we obtain 


‘ ) 
53 (wt+2)(S- VM) +5 (w*+ 302+ 20°)6R | R 
‘ 


+-(§-V6)Vx=0 (23) 

If the perturbation is very small compared to the 
dimensions of the vortex, i.e., if (R-R’)/R=~107', the term 
in 6? in the above equation may be neglected and we 
finally have 


as 2 ee 
dt? R? ” 


Let us assume that the displacement § is proportional 
to a funetion of the form e®! and project equation (24) 
along two orthogonal directions: one corresponding to the 
If the com- 


ay($-M,) ]R+ (§-V6)Va=0 (24) 


vertical and the other to a radial direction. 
ponents of § along these directions are given by 2 and r, 
respectively, the projections of equation (24) along these 


directions are: 


2wa,f OM, , OM, Or/ O68, O86 
ata (" ort? Oe ) aA" a =9:) ’ 
and (25) 
Baal BA. ORE. Or/ O08. O86 
site R (’ or 7 ) ar "Or 25:) 


Where we denotes the absolute rotation (w+). 
lor nontrivial values of 7 and z, the determinant of the 
above system of linear equations vanishes and we obtain 
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(de 29 de B9\ (vg: OM, Ze, OM, 
os “[( Sate a) +( , Er & )] 
2w,, Or ) ( 0M,09 OM, 00 


2w,, Or . 
a: a 0 (2 
( R o)\ dz Or. or 5: ) (<0) 


R Oz 


Examining the order of magnitude of the various terms in 
equation (26), we find that 





5 
2 OM, OM, os 
) —w — — Py). 4 = 8 
R ( oes vr” oz ) 2w-VXV,~10 2 
2 Or Or\ (OM, 08 OM, 00 
R ( ue. det Mer ra) )( dz Or Or >) 
(2w,-Vr)(V8-VX V,) ~10 =a 


Solving for v? and neglecting the last two terms of (27) 
in comparison with the first, we have 


v," Vr-V60 (28) 
The value of this root is determined mainly by the static 
stability 
Or 08 
Oz 02 


q O86 
6 Oz 


(29) 


The second root y,” is obtained by dividing the third 


term of (27) by »,°. Thus, 
2 Or Or\ (OM, 00 OM, 00 
2 hk (4. zt or )( Oz Or Or Oz ) (30 
= Or O9 , On O80 mu) 
Oz Oz Or or 


The above equation will be greatly simplified if it is 
projected on an isentropic surface so that 06/0r vanishes. 
In addition, since the angle between isobaric and isentropic 
surfaces is generally small, (O7/Or) ~0 and 


vq’ =|2wa(V X V).|6 (31) 


Equation (31) is identical with Solberg’s results for the 
polar vortex and states that the square of the frequency 
of the second root is equal to twice the product of the 
absolute rotation of the vortex and the vertical component 
of the absolute vorticity, as measured on an isentropic 
surface. If we now use a relative frame of reference fixed 
with respect to the earth, we have 


Since we are considering a steady vortex, the streamlines 
coincide with the trajectory and R=R, The frequency 
v, becomes imaginary and the motion unstable when 


( ntl) (¢ NI, 0 33) 
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In other words, dynamic instability occurs with normal 
winds provided the absolute vorticity, as measured on an 
isentropic surface, is negative. Alternatively, instability 
occurs with positive absolute vorticity provided the wind, 
as measured on an isentropic surface, is anomalous. 

The above criterion is identical with that obtained by 
van Mieghem [17]. 


visage the possibility of the occurrence of anomalous 


The latter, however, did not en- 
winds and, therefore, equated the condition for the release 
of dynamic instability with the occurrence of negative 
absolute vorticity. 

It should be noted that, according to the theory of 
equations, the second root may be obtained by subtracting 
v, from the negative value of the coefficient of v? in equa- 
then neglect the vertical wind shear 


tion (26). If we 


(O.\M,/0: 


0), we have 


ys? -( 2 +f) (¢+f) (34) 
2 R, ; 

which is similar to equation (32) except that the quantities 
are now measured on a level surface. Equation (34) is 
identical with the criterion for dynamic instability ob- 
tained by Sawyer [15]. 

Physically, the operation of the above criterion may be 
visualized by studying the balance of forces on an anti- 
eyelonic vortex which is initially in gradient balance. 
Let a particle be given an impulse toward lower pressure. 


shown at the 


20 
NAUTICAL MILES 


Vertical cross-section of temperature anomalies from the mean August atmosphere in a direction perpendicular to that of the 


right edge indicate the levels at which data were 


If the flow is unstable, the particle tends to continue in 
the direction of the impulse. There are two cases to 
consider: 


Case I: Instability with Normal Winds. 
normal, there is instability if the pressure gradient along 


If the winds are 


the path of the particle increases at such a rate that the 
speed of the particle remains subgradient. This would 
require a rapid outward increase in anticvelonic rotation 
or a rapid decrease in cyclonic rotation expressed by the 


inequality 


ee 
LR or (I r) |. 


Clase IT: Instability with Anomalous Winds. From figures 
1 and 2, it is seen that if the wind regime is anomalous, 
the speed increases with decreasing pressure gradient and 
trajectory curvature. Therefore, if the trajectory curva- 
ture of the particle remains constant, the pressure gradient 
along the trajectory should decrease at a fast rate; other- 
wise, the speed of the perturbed particle will eventually 
reach a speed greater than that appropriate for gradient 
Actually, 
however, the trajectory curvature decreases as the particle 
Therefore, instability can 


equilibrium and its acceleration will be checked. 


moves toward lower pressure. 
occur with a slower rate of pressure gradient decrease or, 
if the curvature decreases rapidly, instability can oceur 
even with a slow outward increase of the pressure gradient. 
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Figure 11,—Vertical cross-section of deviations of the altitudes of pressure surfaces from the mean August atmosphere (“D’’ values) in 






a direction perpendicular to that of hurricane Daisy, August 25, 1958 (analysis by J. A. Colén). 











Equation (33) states that the motion is unstable so long A comparison of figures 6, 10, and 11 indicates that there 
as the combined effect of the change of wind speed curva- are few solenoids at the level and location where the 





ture is such that anomalous winds occur. We may, therefore, set 


i ew dC. .,4> 7 
> on, 26 a 20 Sand 0 (38) 
la = (VA | s (36) Hh = 8 















. aes ‘ eee aia biieoas C4 2035 constant (39) 
6. A MECHANISM FOR THE DEVELOPMENT 
OF ANOMALOUS WINDS If we approximate the hurricane to a simple circular 





vortex, the circulation at a distance R from the center is 






An insight into the manner in which anomalous winds 






develop may be gained by considering the circulation C=2rR*w (40) 







theorem of Bjerknes which may be written [11]: neal 
dC, 40,45, > >=? sin ¢ (41) 












where @ is the latitude. 






7 Substituting the above values of C and = in (39), we 
Na,- 
obtain 
’ n . ‘ 22,10.) ‘onstar (42) 
OF p V.° dr=absolute circulation around a_ closed R*(w+2,)=constant 4 
~“@ 
curve Q. whence it can be seen that, in a circular vortex, disregard- 





. ing the solenoid term in equation (37) is equivalent to 
; V+ dr=relative circulation around the same curve. assuming conservation of absolute momentum—a concept 
/Q 





which is known to be applicable to upper-level hurricane 





2.=the equatorial projection of the area within the circulation. 






carve ©. According to equation (42), if the vortex has horizontal 
N, the number of pressure-volume solenoids enclosed convergence, the circles contract and w increases. If, on 





by the curve Q. the other hand, the circles expand as a result of horizontal 
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divergence, w decreases. In this case, in order that 


should decrease sufficiently so that 


2V 


R +S - 





4 
R must approach infinity. Since the expansion cannot 
proceed so far, we conclude that the anomalous winds in 
figures 4, 5, and 6 were not produced by a symmetrical 
horizontal divergence of -the air particles which arrived 
at these higher levels by ascent in the core of the storm. 
Another mechanism must be sought to account for their 
occurrence. 

From figure 1 it is clear that a change of regime from 
normal to anomalous winds would be more difficult to 
accomplish if the pressure gradient force is much lower 
than the maximum value /?/4A’,, since this would re- 
quire a big jump in the wind speed. If, however, the 
pressure gradient is nearly equal to the maximum, a 
comparatively small impulse would be sufficient to shift 
an air particle from the lower to the upper branch of the 
) curve of figure 1. 

Let us visualize initial conditions characterized by an 
anticyclonic air stream in gradient equilibrium. In 
figure 12, let A be the equilibrium position of a particle 
in this stream. Now let us visualize that the equilib- 
rium is disturbed by an increase in the pressure gradient 
brought about, for instance, by an interaction between 
tropical and extratropical pressure systems. Two cases 


may be discussed: Case A.-The increase of the pressure 
gradient force is from P to P,, 1.e., 


gradient force is well below the maximum /?/4K%. 


ihe resulting pressure 
The 
wind, having become subgradient, the particle turns 
toward lower pressure and accelerates. Having reached 
the equilibrium speed B, unless there is considerable 
damping, the particle may slightly overshoot this posi- 
tion to By. 
about the equilibrium position until balance is finally 
reached at B. 

Case B. 
from P to P,; i.e., the increased pressure gradient force 
is nearly equal to the maximum value /?/4Aj. In this 
case, the speed of the accelerating particle, by over- 


It then oscillates with decreasing amplitude 


The increase of the pressure gradient force is 


shooting the equilibrium position, may reach the value B; 
which is in excess of the critical speed f/2K). Whereas 
at B, the speed of the particle is supergradient, and at 
Thus, 


by overshooting the critical speed, the particle is hence- 


B, it is subgradient for the same pressure gradient. 


forth constrained to move toward lower pressure and 
sustain further the distribution of 
the wind field along its path is such that the absolute 


acceleration unless 


vorticity is negative; in this case the motion will be 
stabilized in accordance with equation (33). 
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Illustrating the development of anomalous winds 


(see text). 


Figure 12. 


7. EXAMPLES OF SITUATIONS IN WHICH 
ANOMALOUS WINDS PLAY A ROLE 


In a previous paper [1], the author suggested from 
theoretical considerations that the occurrence of anoma- 
lous winds in the manner described above provides a 
De- 


tails concerning how this is brought about and some 


mechanism for triggering hurricane development. 


observational evidence to this effect are given in a sepa- 
rate paper [2] and are beyond the scope of the present 
study. Here we would suggest that, in view of their essen- 
tially unstable nature, anomalous winds must play an im- 


portant role not only in hurricane formation but also in 


other types of atmospheric development, especially those 
which occur downstream from a pressure ridge which 
has undergone marked intensification, or those accompany- 
ing the superposition of strong westerly winds on the 
crests of the long waves. 

A case in point is that discussed by Bjerknes [4] linking 
the deepening of a wave trough with a strong intensifica- 
tion of the ridge upstream from it. Specifically, Bjerknes 
suggests that if the ridge intensifies so that the curvature 
of the contours exceeds a critical value 


44 


where V, denotes the geostrophic wind, the air particles 
cannot follow the contours and must perforce flow to 





inds 


Novi MBER 1961 


Figure 13.—-Illustrating the cyclogenetic effect of anomalous winds. 
The dashed curve represents the probable trajectory of an air 
particle after its speed changes to the anomalous regime at the 
point A, 


ward lower pressure and accelerate. The acceleration 
results in supergradient winds which then curve back to 
higher pressure, the net result being a deepening of the 


trough downstream from the ridge. 


Now that the likelihood of the occurrence of anomalous 
winds in situations analogous to that described above has 
been demonstrated, the deepening of the trough may be 
looked at in a new light. Let us suppose that the wind 
becomes anomalous as a result of the strengthening of 
the pressure gradient to a value near the maximum pos- 
sible for anticyclonic flow, and let us follow the trajectory 
of a particle situated at the point A in figure 13, where the 
wind changes regime and becomes subgradient in the 
As a result, 
the air particle moves to lower pressure and accelerates. 


manner discussed in the previous section. 


In doing so, its curvature decreases, as can be seen from 
the figure. Unless the pressure gradient along the tra- 
jectory increases rapidly, the progressively decreasing 
curvature requires an increasingly stronger speed for 
balanced motion (see fig. 2) and the particle, therefore, 
However, if the 
vanishes so that the particle follows a straight trajectory, 


continues to accelerate. curvature 
the wind instantly becomes supergradient and the particle 
must curve back to higher pressure. In the mean, the 
particle must, therefore, follow a slightly anticyclonic 
trajectory and is constrained by the geometry of the 
pressure distribution to curve back to higher pressure, 
as shown by the dashed curve in figure 13. 


In this connection it is of interest to mention that 
Angell [3] has noted that anomalous flow toward higher 
pressure can occur continuously for several hours. Thus, 
by stipulating the occurrence of anomalous winds, we 
arrive at the same ultimate deepening envisaged by 
Bjerknes. 
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8. CONCLUSION 


Although the observational evidence presented here 
is for the occurrence of anomalous winds over small regions 
of the atmosphere, it is significant to note: 

A. that their occurrence may be dependent on large- 
scale processes, such as those involving the intensification 
of pressure ridges or the latitudinal shift of the westerlies; 
and 

B. that their effect extends beyond the area where they 
occur. 

The fact that these winds can occur only in anticyclonic 
motion again demonstrates the vital regulatory role of 


the high pressure regions in determining the state and 


motion of the atmosphere. 
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THE WEATHER AND CIRCULATION OF AUGUST 1961 


Record Heat from the Northern Plains to the Pacific Coast 
JAMES F. ANDREWS 


Extended Forecast Branch, U.S. Weather Bureau, Washington, D.C, 


1. HIGHLIGHTS 


August 1961 continued the succession of hot dry months 
in the northwestern United States! with average and 
daily temperatures reaching record levels in many areas. 
Severe drought persisted in the Northern Plains and by 
the end of the month conditions were comparable with 
the most severe drought of other years. 

The Midwest and South were unseasonably cool, con- 


tinuing a pattern that began in April of this year. The 


most persistent coolness occurred in the western Gulf 


States where mean temperatures for the month were at 
or near record levels. 

This sharp contrast in weather regimes was also char- 
acteristic of the summer season and was related to an 
abnormal mid-tropospheric circulation. 


2. MONTHLY MEAN CIRCULATION 


The average circulation pattern at 700 mb. for August 
1961 (fig. 1) featured an unusually strong ridge over west- 
ern North America. The height anomaly center asso- 
ciated with this ridge was the highest ever observed in this 
area during the Extended Forecast Branch period of record 
1933). 
deeper than normal trough over the southeastern United 
States. 

In the Pacific and North America the pattern of 700-mb. 
height and its anomaly exhibited a meridional character, 
with a deep mid-Pacifie trough helping to support the 
strong ridge in the climatologically favored area over 
western North America (fig. 1). Elsewhere in the hemi- 
sphere, from the Atlantic across Asia, the circulation was 
markedly zonal. 

The middle-latitude 
normal and displaced northward during August over the 
Greatest 


(since A companion feature of this ridge was a 


westerlies were stronger than 
western portion of the hemisphere (figs. 2, 3). 
displacement was over North America where a band of 
wind speeds as much as 5 m.p.s. above normal was 
observed across central Canada (fig. 3B). These fast 
westerlies were related to the deep cyclonic center in 
northeastern Canada and strong ridge to the west (fig. 
i). In the United States, however, weaker than normal 
westerlies predominated. The contrasting wind speed 


Reference to the United States in this paper will not include Alaska or Hawaii. 


anomaly pattern in North America (fig. 3B) is also in- 
herent in the 700-mb. height anomaly pattern (fig. 1). 
Across Canada the anomalous flow was westerly, while in 
the United States easterly anomalous flow prevailed. 

The westerlies over the eastern portion of the hemi- 
sphere were displaced south of their normal August 
position, primarily as a result of the blocking High and 
area of 700-mb. height anomaly centered near Novaya 
Zemlya (figs. 1, 3). Fast westerly flow prevailed across 
the Atlantic, Europe, and into Asia, with greatest wind 
speed departures, as much as 6 m.p.s., in central Europe 
(fig. 3B). 


3. AVERAGE MONTHLY WEATHER 
RECORD HEAT AND DRYNESS IN THE NORTHWEST 


The northwestern United States, from the northern 
Great Plains to the Pacific coast, experienced one of the 
hottest and driest months ever recorded. Temperatures 
for August averaged from 4° to more than 6° F. above 
normal in all but coastal areas (fig. 4), while less than half 
the normal amount of precipitation fell in the Northern 
Plains (fig. 5). In many areas stations with long-period 
records reported that August 1961 was the warmest month 
of record. <A representative list of these stations is shown 
in table 1. 

The hottest weather occurred quite generally during 
the first week when a strong ridge dominated western 
North America. Extreme heat from the Cascade Moun- 
tains to the northern Great Plains sent temperatures well 


over 100° F., with many stations establishing new absolute 


Monthly mean high temperature records established in 


August 1961 


TABLE 1. 


Temperature (°F.) 
Station ‘i 
Year records 
began 


Departure 
from normal 


Average 
monthly 


Devils Lake, N. Dak 72.6 5 1904 
Williston, N. Dak ‘ 7 1878 
Rapid City, 8. Dak 77 : 1887 
Billings, Mont 76. 2 5 1804 
Glasgow, Mont 5 1843 
Great Falls, Mont 5.3 Before 1900 
Havre, Mont 7 t 1879 
Sheridan, Wyo 7 5 1907 
Seattle, Wash 9. 1892 
Portland, Oreg 72.2 : 1902 
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Figure 1.—Mean 700-mb. contours (solid) in tens of feet, and departures from normal (dotted) in 50-ft. intervals, for August 1961. 
Important circulation features were the unusually strong ridge over western North America and deep trough over southeastern United 


States. 


TABLE 2.—Mazximum temperatures 


during August 


1961 


which 


equaled or exceeded previous highest recorded for any month. 


Station 


Helena, Mont 
Kalispell, Mont 
Missoula, Mont 
Havre, Mont 
Boise, Idahe 
Spokane, Wash 
Pendleton, Oreg 


*Equaled previous record 


Temperature 
F 


*103 
105 
105 
111 
110 

*108 


113 


Date 


maximum temperature records as shown in table 2. In 
addition to the stations listed, a record high August 
temperature of 115° F. was reported at Lewiston, Idaho 
on the 4th. Nights were also unusually hot in the north- 
ern Plains where minima were occasionally in the mid-70’s. 

August continued the series of dry months that has 
left soil moisture reserves at a critically low level. At 
Williston, N. Dak., only 10 percent of the normal moist ure 
fell, the driest August since observations began in 1878. 
Much of the precipitation in the Northwest was in the 
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Ficure 2.— Mean 700-mb. zonal wind speed profiles in the western 
portion of the Northern Hemisphere for August 1961 (solid) and 
August The middle-latitude westerlies were 
stronger than normal and displaced northward. 


normal (dashed). 


form of widely scattered light showers which failed to add 
to soil moisture reserves. 

Drought conditions in the Northwest were related to 
the abnormally strong ridge and area of positive 700-mb. 
height anomaly over western North America (fig. 1). 
Advection of voriicity from the deep mid-Pacific trough 
strengthened this ridge, normally found over western 
North America during the summer months. The wester- 
lies, displaced some 10° north of their normal position 
over western Canada (fig. 3A), diverted cool, moist 
Pacifie air masses well to the north of the drought area. 
At the end of the month, however, progression of the 


long-wave pattern began, thus bringing an abrupt end to 
the heat wave in the Northwest along with the first sub- 
stantial rains of the month to Washington and northern 


Oregon. 
PERSISTENT COOLNESS IN THE MIDWEST AND SOUTH 


This was the fifth consecutive month that temperatures 
have averaged below normal in the area from the central 
Plains to the Appalachians and southward to the lower 
Mississippi Valley (fig. 1 in [1] and fig. 4). Greatest de- 
partures during August, as much as —4° F., were observed 
in eastern portions of the central and southern Plains. 
This was the coolest August of record at Shreveport, La., 
Austin and Midland, Tex. Daily temperatures 
averaged below normal on nearly every day of the month 
in the west Gulf States with some cities, such as Jackson, 
Miss. and Shreveport, La., never exceeding normal. 

These cool conditions were related to the trough over 
the Southeast and its associated large area of negative 
700-1nb. height anomaly (fig. 1). Furthermore, this cireu- 
lation favored the unusually far southward displacement 
of cold fronts with their attendant cool air masses. 


and 
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Figure 3.—(A) Mean 700-mb. isotachs, and (B) departures from 
monthly normal wind speeds, both in meters per second, for 
August 1961. Solid arrows in (A) indicate principal axes of 
maximum wind speeds, and dashed arrows their normal August 
positions. The westerlies over North America were displaced 
north of normal. 


The distribution of precipitation bore a nearly typical 
relationship to the 700-mb. circulation. Less than normal 
amounts fell generally west of the mean trough while 
greater than normal amounts fell east of the trough in 
stronger than normal southerly flow (figs. 1, 5). Following 
two months of excessive precipitation, central Texas was 
unusually dry with the Dallas-San Antonio areas receiving 
less than 10 percent of their normal amounts. Frequent 
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Figure 4.--Departure of average surface temperature from normal 
F.) for 1961. 
prevailed in the Midwest. 


A strong temperature anomaly gradient 
(From [8] 


August 


showers in the South Atlantic States resulted in many 
rainy days (0.01 inch or more) with Charleston, 8.C., 
having 19 such days, the greatest number in August since 
1898. 

WEATHER ELSEWHERE IN THE NATION 

A warm wet weather regime prevailed in the Southwest 
during August. Temperature departures were rather 
small, +2° F. or less in most areas (fig. 4), while precip- 
itation amounts were mostly greater than normal. 

The Southwest moist tongue was unusually active with 
shower and thunderstorm activity very frequent. Ely, 
Nev., had 24 days with the greatest 
number of days with thunderstorms ever recorded there 
Some of this activity spread 


thunderstorms, 


in any month in any year. 
into California where precipitation in August is very 
uncommon. The warin showery weather in the South- 
west Was associated with above normal mid-tropospheric 
heights and easterivy anomalous flow (fig. 1). 

The Northeast normal with 
temperature departures mostly less than 2° F. (fig. 4), 
while precipitation was generally near normal. In New 


England, the first 20 davs of the month were very dry, 


was also warmer than 


the result of westerly anomalous flow in mean trough 
During the latter part of the 
As a result 


conditions over the area. 
month this trough retrograded and weakened. 
south and southeasterly anomalous flow became estab- 
lished, accompanied by frequent periods of precipitation. 
Most of the month’s total rainfall fell at this time. 


4. TROPICAL STORM ACTIVITY 


During any August one or more tropical storms can be 
expected to develop in the Caribbean or North Atlantic [2], 
but none occurred in August 1961. It is generally con- 
sidered that tropical storm development is related to the 
strength and position of the middle-latitude westerlies [3]. 


During this August the average circulation appeared 
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(inches) for August 1961. The 


(From [8].) 


Total 
Northern Plains were very dry. 


FiGcure 5. precipitation 


quite favorable for generation of such storms. The 700- 
mb. height anomaly pattern agreed quite closely with the 
findings of Ballenzweig [4] for maximum storm occurrence 
Furthermore, the associated upper-level westerlies were 
displaced north of normal (fig. 3A). It is evident that 
some of the basic elements necessary for tropical storm 
development were missing. 

Tropical storm activity was confined primarily to the 
western Pacific, where four storms developed, three of 
which reached typhoon intensity. 

Typhoon June developed near 12° N., 133° E. 
Ist. This storm moved steadily northwestward, losing 


on the 


much of its strength as it passed over Taiwan on the 6th, 
finally dissipating after moving into China on the 8th. 
Typhoon Kathy first appeared near Iwo Jima on the 15th. 
It also moved northwestward, striking Kyushu, Japan, on 
the 17th before filling on the 18th. The largest and most 
intense storm was typhoon Lorna. 
16° N., 129° E. on the 20th, Lorna grew rapidly to typhoon 
intensity. Moving northwestward Lorna’s “eye’’ passed 
over southern Taiwan on the 24th, reached the China 
coast on the 25th, and dissipated inland on the 26th. 
Tropical storm Marie developed near 23° N., 153° E. 
on the 29th. This storm moved slowly northward, turned 
on a westerly course and dissipated early in September, 


First appearing near 


never reaching typhoon strength nor threatening any 


land areas. 


5. THE SUMMER IN REVIEW 


Summer of 1961 in the United States was characterized 
by contrasting regimes of weather associated with an 
abnormal mid-tropospheric circulation. In figure 6 are 


shown (A) the average circulation at 700 mb., (B) the 


departure of average surface temperature from normal, 
and (C) the total precipitation, for the months of June, 
July, and August. 
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Fiaure 6.—(A) Mean 700-mb. contours (solid) in tens of feet and 
departure from normal (dotted) in 30-ft. intervals with centers 
labeled in tens of feet, (B) Departure of average temperature from 
normal (°F.), and (C) total precipitation (inches). All charts 
for summer (June through August) 1961. An abnormally strong 
ridge dominated the circulation over North America bringing 
severe drought to the Northern Plains. (B and C from [8}.) 


Unseasonable warmth prevailed throughout most of 
the western half of the Nation with greatest temperature 
departures, 4° F. or more, in the Northwest (fig. 6B). 
This was the warmest summer ever recorded in Walla 
Kalispell, 
Mont., with a record dating to 1899, also had its warmest 


Walla, Wash. since observations began in 1886. 
summer. Near record warmth was observed elsewhere 
in the Northwest. 

Extreme dryness also prevailed in the area of greatest 
At Billings, Mont. this was the driest summer 
Devils Lake, N. 


second driest summer. 


warmth. 
on record, while Dak. experienced its 
This hot dry weather regime was accompanied by 
the Northern Plains, conditions 
severe in northwestern North Dakota and 
An index to measure drought severity 
ry’ ° 

The index 


increasing drought in 
being most 
eastern Montana. 
has recently been developed by Palmer [5]. 
values depend on the temperature and precipitation during 
each month, the weather during preceding months, the 
climate of the place being studied, and the duration of the 
abnormally dry weather. At the end of May 1961 the 
drought index for northwestern North Dakota was —0.60 
incipient drought stage). By the end of August the 
had 4.20, comparable to the extreme 
drought of 1934 when the drought index was —4.17 at 


index reached 


—_—. 


* A value of —4,00 is considered to be the threshold between severe and extreme drought, 


























the end of August. Crop correspondents of the Depart - 
ment of Agriculture, in their report of September 11, 
1961 


existed in eastern Montana and the western Dakotas 


[6], considered that extreme drought conditions 

The tendency for summer drought to follow a warm 
dry spring in the Great Plains has been shown by Namias 
(7). 


was warm and generally dry in the Northern Plains [8]. 


This vear appears to be no exception since spring 


Furthermore, these conditions have existed since summer 
1960, as is substantiated by the Palmer index which had 
been at the incipient to mild drought stage since that 
time. 

In marked contrast with the West, most of the eastern 
half of the Nation enjoyed a cooler than normal summer. 
Greatest temperature departures (as much as 4° F. in 
some areas) and near record coolness were observed in 
the South (fig. 6B). It 
reported at Baton Rouge, La. in a record dating to 1893, 


was the coolest summer ever 
while at Abilene, Tex., this was only the second summer 
since 1886 that the temperature has failed to reach 100° F. 

Precipitation averaged somewhat greater than normal 
in much of the East and South, with a large area in eastern 
and central Texas receiving more than two to three times 
the normal amount (fig. 6C). Most of the latter pre- 
cipitation fell during June [9] and July [1]. 

Comparison between the mean 700-mb. circulation pat- 
tern for August 1961 (fig. 1) and that for summer of this 





500 


vear (fig. 6A) shows them to be remarkably similar in the 
United States. It follows that persistence of the circu- 
lation and weather anomalies was also very high. It is 
well known that surface temperature anomalies and mid- 
tropospheric height anomalies are most directly related 
in summer. Note the close “fit”? during the summer of 
1961 (figs. 6A, B). The abnormally strong ridge over 
western North America, with its associated subsidence 
and northward displaced westerlies tended to inhibit 
precipitation in the Northwest, with less than half the 
normal seasonal amounts falling over the Northern Plains 
(fig. 6C). An examination of the file of 700-mb. mean 
summer charts in the Extended Forecast Branch discloses 
that the positive height anomaly center over western 
North America during this summer was the strongest and 
most extensive positive anomaly ever observed in this 
area. Retrogression and deepening of the trough nor- 
mally found along the east Coast in summer was effected 
primarily by the unusually strong ridge to the west. 
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Weather Note 
EXCESSIVE RAINFALL IN THE MONTGOMERY, ALA. AREA ON AUGUST 31, 1961 


ARTHUR R. LONG 


Weather Bureau Airport Station, Montgomery, Ala. 


{|Manuscript received September 8, 1961] 


A band of heavy thundershowers occurred in the Mont- 
gomery, Ala. area on the night of August 30-31, 1961. 

At Montgomery the 65-year maximum 1-hour precipita- 
tion record of 3.81 in., which occurred on July 20-21, 1899, 
was broken when 4.16 in. fell between 3:10 a.m. and 4:10 


am. CST, on August 31, 1961. For the 2-hour period 


3:10 a.m. to 5:10 a.m., 4.51 in. fell, breaking the former 


2-hour record of 4.38 in. set on November 9, 1927. 


Several cooperative stations received heavy rainfall also, 


most of which fell between midnight and 7 a.m. on August 


31. Some of the heavier amounts were: Mathews, 6.05 


in.; Troy, 6.75; Brundidge, 4.13; Elba, 6.58; and Geneva, 
5.64 in. At Milstead, Ala., a total of 11.50 in. fell during 
the 24-hour period ending at 11 a.m. on the 31st. 

This excessively heavy rainfall caused considerable dam- 


age by flooding streets and residences in Montgomery and 


washouts on highways and one railroad in the surrounding 
area. Numerous highways were closed for periods of one 


to three days. 
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